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Abstract 

Recent measurements of the Gottfried Sum Rule have focused attention on the 
possibihty of substantial breaking of flavor symmetry in sea quark distributions of 
the proton. This has been confirmed by by pp and pD Drell-Yan processes measured 
at FNAL. The theoretical models used to infer flavor symmetry breaking rely on the 
assumption that parton distributions are charge symmetric; it is conceivable that 
current tests of flavor symmetry could be affected by substantial charge symmetry 
violation. Since all phenomenological parton distributions assume the validity of 
charge symmetry, in this paper we examine the possibility that charge symmetry 
is violated [CSV]. We first list definitions for structure functions which do not 
make the usual assumption that parton distributions obey charge symmetry. We 
then give some simple model estimates of CSV for both valence and sea quark 
distributions. Next, we list a set of relations which must hold if charge symmetry is 
valid, and we review the current experimental limits on charge symmetry violation 
in parton distributions. We then propose a series of possible experimental tests of 
charge symmetry. The proposed experiments could either detect charge symmetry 
violation in parton distributions, or they could provide more stringent upper limits 
on CSV. We discuss CSV contributions to sum rules, and we propose new sum rules 
which could differentiate between flavor symmetry, and charge symmetry, violation 
in nuclear systems. 
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1 Introduction 



It has long been recognized that the strong interaction respects charge symmetry to a 
high degree. We review the definition of charge symmetry, which is sometimes confused 
with isospin symmetry. For details we refer the reader to comprehensive reviews on charge 
symmetry by Miller, Nefkens and Slaus and Henley and Miller 0. The assumption 
of isospin independence of hadronic forces requires that the Hamiltonian of the system 
commutes with the isospin operator T, i.e. 



[H,T] 



(i: 



Whereas isospin symmetry requires invariance of the Hamiltonian with respect to all 
rotations in isospin space, charge symmetry requires invariance only with respect to ro- 
tations of 180° about the T2 axis, where the charge corresponds to the third axis in 
isospin space. Consequently, isospin symmetry necessarily implies the validity of charge 
symmetry; however, the converse is not necessarily true. As charge symmetry is a more 
restricted symmetry than isospin symmetry, it is generally conserved in strong interac- 
tions to a greater degree than isospin symmetry. Thus, while in many nuclear reactions 
isospin symmetry is violated at the few percent level, in most cases charge symmetry is 
obeyed to better than one percent. 

For a system of A particles, the charge symmetry operator can be written as 

Pc. = exp(mT2) = Hf^iexp[mr2(0] . (2) 

The operation of charge symmetry maps up quarks to down, and protons to neutrons. 
Specifically, under charge symmetry 

e*'"^^ : u ^ d,p n 

: d —u, n — > —p (3) 

At the quark level, charge symmetry implies the invariance of a system under the 
interchange of up and down quarks. The proton and neutron each contain three valence 
quarks, plus a "sea" of quark-antiquark pairs. Coulomb effects aside, the "proton" is 
converted to a "neutron" by interchanging up and down quarks in the two nucleons. At 
the level of parton distributions, charge symmetry implies the relations 

dP{x,Q') = m"(x,Q2) (4) 

Charge symmetry is broken by electromagnetic interactions, but these should play a 
minor role at high energies. We can get an estimate for the magnitude of charge symmetry 
violation [CSV] at the parton level; we would naively expect parton CSV to be of the order 
of the up-down current quark mass difference divided by some average mass expectation 
value of the strong Hamiltonian, or (m^ — m„)/ < M >, where < M > has a value of 
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roughly 0.5-1.0 GeV. This would naturally put CSV effects at a level of 1% or smaller. 
Note that we expect charge symmetry to be valid at this level, despite the fact that the 
current quark masses themselves, i.e. ~ 4 MeV, ~ 7 MeV, differ by 50%! However, 
our understanding is that dynamical chiral symmetry breaking and/or confinement masks 
this very large "primordial" violation of charge symmetry, and observable quantities are 
expected to respect charge symmetry to roughly one percent. 

In nuclear physics, charge symmetry involves the interchange of protons and neutrons 
in a system. At low energies, charge symmetry appears to be generally valid at the level 
of 1% or better in nuclear systems, although there are some notable exceptions to this 
rule of thumb The proton and neutron masses are equal to about 0.1%; the binding 
energies of tritium and ^He are equal to 1%, after Coulomb corrections. We can compare 
energy levels in "mirror" nuclei (nuclei related to one another by Pes), and generally find 
agreement to better than 1%, after correcting for electromagnetic interactions. 

From our experience with charge symmetry in nuclear systems, and because of the 
order of magnitude estimates of CSV in parton systems, charge symmetry has been 
universally assumed in quark/parton phenomenology. With this assumption, one reduces 
the number of independent quark distribution functions by a factor of two. One simply 
defines all quark distribution functions in the neutron to be equal to the corresponding 
functions in the proton, while interchanging up and down quarks in the process. The 
assumption of charge symmetry is sufficiently ingrained in quark/parton phenomenology 
that its validity is a necessary condition for many relations between structure functions. 
Thus, it is not apparent to many physicists that several sum rules or structure function 
equalities may be valid only to the extent that charge symmetry is exact. 

Recently, much attention has been focused on the apparent violation of what is called 
SU(2) flavor symmetry in the nucleonQ. The measurements of the Gottfried sum rule 
reported by the New Muon Collaboration (NMC) [Q sparked a great deal of interest in the 
sea-quark flavor distributions of nucleons [^HTTH. The "natural" explanation of the NMC 
results is that "SU(2) flavor symmetry" is broken in the proton sea quark distributions 
(i.e. u^{x) 7^ d^{x)). This has been widely cited, and several theoretical investigations 
have been carried out to investigate the possible origin of this flavor symmetry violation. 
At the time of the NMC measurements, another possible explanation for this data was 
that the Gottfried sum rule was obeyed, and that the apparent experimental violation 
resulted from significant contributions to the sum rule at extremely small values of x |T^ . 

Since the NMC experiment suggested large flavor symmetry violating antiquark dis- 
tributions, this led people towards experiments which had the possibility of "direct" 
observation of flavor symmetry violation [FSV] in the proton sea. Ellis and Stirling [|1^ 



pointed out that this information could be obtained by studying proton-induced Drell- 
Yan cross sections with both proton and neutron [i.e., deuteron) targets. Two subsequent 



experiments have been carried out, by the NA51 group at CERN |]14|, and the E866 ex- 



periment at FNAL [ITH] . Experimental results from these collaborations also seem to show 



^We adopt this terminology, which is widespread, despite the fact that there is no underlying SU(2) 
symmetry which is broken. The term refers to the fact that proton sea quark distributions are not equal, 
i.e. dP and are not equal at all Bjorken x. 
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a large flavor symmetry violation [FSV] in the proton sea quark distributions. We will 
discuss these results in detail in Sect. 4 of this review. 

However, the conclusion that these three experiments demonstrate large FSV effects, 
as well as the magnitude of the FSV effects extracted, relies on the implicit assumption 
of charge symmetry. It has been pointed out [16, l^] that all three experimental results 
(NMC, NA51 and E866) could in principle be explained even if flavor symmetry were 
conserved, if we assume charge symmetry violation [CSV] in the nucleon sea. As we will 
show, the CSV terms necessary to account for the results of these experiments would 
be surprisingly large - much larger than theoretical models would predict. If charge 
symmetry were violated to this degree in parton distributions, it would be amazing that 
at low energies it would be nearly an exact symmetry. 

However, the history of our understanding of nucleon structure has involved a series 
of similar "surprises" from experiment: e.g., the significant contribution of antiquarks 
to nucleon structure functions at small x; the large fraction of the nucleon's momen- 
tum carried by glue; the persistence of substantial spin effects at high energies despite 
perturbative QCD [pQCD] predictions that single-spin asymmetries should vanish; the 
"spin crisis," which suggests that a surprisingly small fraction of the proton's spin may 
be carried by valence quarks; and the behavior of quark distributions at very small x. So, 
despite the strong indirect evidence from low-energy physics, and straightforward pQCD 
arguments which suggest that charge symmetry should be valid to about the 1% level in 
parton distributions, we urge the reader to keep an open mind on this question. 

In this paper, we will provide a comprehensive review of the following question: how 
valid is the assumption of charge symmetry for parton distributions? First, we will 
redefine the nucleon structure functions in terms of quark/parton distributions, without 
assuming charge symmetry. Next, we will show how relations between structure functions 
become modified when we allow CSV terms. We then calculate the CSV contributions to 
various observables. We examine the current experimental evidence for charge symmetry. 
As we will show, all experiments to date are consistent with parton charge symmetry. 
However, in some regions present experimental upper limits on parton charge symmetry 
violation are rather weak. On the other hand, new experimental neutrino deep inelastic 
scattering data, when taken together with high energy muon scattering, can provide 
rather strong constraints on parton CSV, at least for a certain range of Bjorken x. 

We will also present some simple model estimates of charge symmetry violation in 
both valence quark and sea quark distributions. For the "majority" valence quark distri- 
butions, {i.e., — we predict very small CSV amplitudes, no larger than 1%. 
However, our model calculation of charge symmetry violation in the "minority" nucleon 
valence quark distributions {d^{x) —u^{x)) [jl8|,|l9| suggests surprisingly large CSV terms. 
We discuss several experiments which could detect CSV in parton distributions, or which 
could improve the current upper limits on quark CSV. 

The structure of our paper is as follows. In Sect. 2 we review the general expressions 
between cross sections and structure functions for deep inelastic scattering processes. We 
write down the most general form of the structure functions, without assuming charge 
symmetry. In Sect. 3, we give derivations, from simple models, of charge symmetry 
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breaking for both valence quarks and sea quarks. We show the magnitude and sign of 
the expected CSV terms in these models. We also review relations between structure 
functions which hold if charge symmetry is valid. 

In Sect. 4 we review those experiments which currently place the best upper limits on 
CSV in parton distributions. Because of the current interest in flavor symmetry violation 
in the proton sea, and because current "tests" of FSV in fact are testing a combination 
of FSV and CSV, we review at length the recent Drell-Yan measurements which are 
presented as evidence for FSV. We review the constraints which recent experiments place 
on CSV and FSV in antiquark parton distributions. Preliminary results from the E866 
Drell-Yan experiment suggests that they can measure the relative magnitude of d and 
u in the proton, over a fairly wide kinematic region. In this same general region, we 
also have data from the NMC measurement of F2 structure functions in protons and 
neutrons, using high energy muon beams. In addition, we have the structure function 
measured by the CCFR group, from charge changing weak interactions induced 
by neutrinos and antineutrinos on iron. All three experiments obtain measurements at 
similar values of Q"^ and x. 

In Sect. 5, we propose experiments which could in principle reveal charge symmetry 
violation in the valence quark distributions (these would also differentiate between FSV 
and CSV effects). In Sect. 6, we review QCD sum rules. We show how these are modified 
if we include sea quark CSV contributions. We review the best known unpolarized sum 
rules, the Gottfried, Adler and Gross-Llewellyn Smith sum rule. In Sect. 7, we show that 
by defining two new sum rules, it would be possible to measure separately CSV, and FSV, 
contributions to sea quark distributions. We call these the "charge symmetry" and "flavor 
symmetry" sum rules, respectively. We also review the status of existing experiments to 
determine current upper limits on sea quark CSV via the charge symmetry sum rule. In 
Sect. 8 we present our conclusions. 

2 Relations Between High Energy Cross Sections and 
Parton Distributions 



2.1 General form of high energy cross sections 



We can write the cross sections for deep inelastic scattering in terms of a set of struc- 
ture functions, which depend on the relativistic kinematics of the reaction. Through 
the quark/parton model, these structure functions can in turn be written in terms of 



quark/parton distributions |20]. For example, the most general form of the cross section 
for charged current interactions initiated by charged leptons on nucleons has the form 



dx dy 



a 



,2sin2 OwM^^ 



xy'rrM' 



6 




Figure 1: Schematic picture of deep inelastic scattering of charged leptons from a nucleon. 

a) Charged- current weak interactions. An intermediate W is absorbed on the nucleon. 

b) Neutral-current electroweak interactions. 



+ (i - 1/ - ^) Fr {X, Q') Tiy- yV2)xFr {X, Q') 



(5) 



This process is shown schematically in Fig. ||a. It involves a charged virtual of 
momentum q being interchanged between the lepton/neutrino vertex, and the hadronic 
vertex. The relativistic invariants in Eq. ^ are = — g^, the square of the four momen- 
tum transfer for the reaction, x and y. For four momentum k {p) for the initial state 
lepton (nucleon), we have the relations 

y 



X 



2p ■ q 



p ■ k 



(6) 



In Eq. ^, Mw is the mass of the charged weak vector boson, and 9w is the Weinberg 
angle. 

Similarly, the cross section for charged current interactions initiated by neutrinos or 
antineutrinos on nucleons has the form 



cc 



a 



dx dy 



ITS 



2 sin^ OwM? 



w , 



xy'Fr^ix,Q'' 



+ [i-y- 



Fr{x,Q')±{y-y'mxFr{x,Q 



(7) 



This process is obtained by interchanging the initial and final state leptons in Fig. |I]a. 
Neutral current (NC) reactions initiated by neutrinos or antineutrinos have the form 

/ \ 2 / , ^9 \ 2 



d'a 



NC 



a 



dx dy 



ns 



Ml 



2 sin^ Ow cos^ OwMl 



Ml 



xy'Ff{x,Q' 
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+ {l-y- ^) Fiix, Q') ±{y- yV2)xFiix, Q^) 

Finally, the cross section for scattering of a left (L) or right (R) handed charged lepton 
in NC reactions has the form 



dx dy 



Allans 



{\xy'F2{x,Q') + {l-y)F:j{x,Q') 



ve ± ai 



(g2 + M|) 2 sin cos^ 



± {y-yy2)xFS^{x,Q 



w 



+ 



xy^Ff{x,Q^) + {l-y)Ff{x,Q' 



(g2 + M|)y 2 sin cosily 
+ (1 - y)F^{x, Q^) ±{y- y^/2)xFi{x, Q') 



xy'Ff{x,Q^ 



(9) 

This process is shown schematically in Fig. |I|b. Either a photon or boson can be 
exchanged in this process. 
In Eq. ^, we have 

-1 + 4sin2 9w 
4 sin 6w cos 6w 

ciiij„ — dip 



4 sin 9w cos 9w 



(10) 



Eq. 1^ describes the deep inelastic scattering for an L (K) handed charged lepton from a 
nucleon. For momentum transfers which are sufficiently small (relative to M|), we can 
neglect the contribution from Z bosons, in which case the scattering is a function only 
of the two electromagnetic structure functions, and F2, respectively. 

2.2 Structure functions in terms of quark/parton distributions 

The form of the proton structure functions, obtained from deep inelastic scattering of an 
electron or muon, can be written in terms of interaction of the charged leptons and quarks 
with the virtual photon 7 Here we assume we are at sufficiently low energies that 
we can neglect contributions from Z^ in electroweak processes. From Eq. |^ we see that 
the resulting cross section can be written in terms of two structure functions, Fi and 
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Furthermore, we work in an energy regime where both and the energy transfer are 
very large, while x remains finite, so that scaling is valid, i.e. the structure functions (to 
first approximation) depend only on x and not on Q^. The resulting structure function 
Fi^{x) can be written in terms of the parton distributions as 



F^^ix, Q^) = -(- yix) + u'ix) + c^ix) + c^ix)] 

2 v y 

+ I \dP{x) + dF{x) + s^{x) + s^{x) 



9 



(11) 
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Figure 2: Coupling of a charged lepton to quarks through exchange of a virtual photon. 



This process is shown schematically in Fig. In Eq. |TT|, we assume we can neglect any 
contribution from bottom or top quarks in the proton. The virtual photon couples to the 
squared charge of the struck quarks. To obtain the corresponding Fi structure function 
for the neutron, we simply change the superscript p ^ n everywhere in Eq. |ll[ 

In Eq. |Tl| (and in most subsequent equations), we have neglected the dependence 
of the parton distributions on the scale at which they are evaluated. As is well known 
2^,^, there is an uncertainty in the parton distributions with respect to the scale /i^ 



at which they are evaluated. Once we calculate the parton distributions qk{x,^'^) and 
gluon distributions at some starting scale, we can evolve the parton distributions to some 
higher through the QCD evolution equations of Dokshitzer, Gribov, Lipatov, Altarelli 
and Parisi For convenience, we will generally omit this scale in equations involving 



parton distributions. 

In the lowest order quark/parton model, the structure function F2 is related to the 
structure function by the Callan-Gross relation |^ 

F2^(x,g2)=2xF7(x,Q2) (12) 

As is well known, the Callan-Gross relation is valid if the virtual photon which initiates 
this process is completely transverse. The more general relation between the two structure 
functions is 

In Eq. 0, i? = (Tl/ctt is the ratio of the cross section for longitudinally to transversely 
polarized photons. An analogous relation will hold for the weak structure functions Ff^ . 
An empirical relation fit to the world's available data on R has been made by Whitlow 
et al. 1^^. The formula is 



ln(QV0.04) g2 g4 + o.09 
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Figure 3: a) Quark contributions to the charged-current reactions induced by neutrinos. 
The virtual is absorbed by negatively charged quarks, b) Antiquark contributions 
to neutrino induced charged-current reactions. 



e 



1 + 



(14) 



The coefficients in Eq. 14 can be found in Ref. |24|. This fit covers the region accessible 
at that time, i.e. x > 0.1 and < 125 GeV^. 

Charged current neutrino scattering on hadrons is mediated by emission of the weak 
vector boson W"^ by the leptons and subsequent absorption of the W"^ on the proton or 
neutron. Thus the structure function Fi corresponding to charge- changing interactions 
of neutrinos on protons can be written in terms of the quark distribution functions as 



+ cP(a)|Kfer^(x6 -x) + c*(x)[|Kd|' + \Vcs\% 



(15) 



In Fig. ^ we show the coupling of the virtual W'^ to quarks; the coupling is to quarks 
with negative charge. In Fig. ^ we show the coupling to antiquarks. 

Similarly, the structure function Fi corresponding to charge- changing reactions for 
antineutrinos on protons can be written as 



Fr ^(x) = MP(x)[iKdr + iK.ii+Mne6)iK;,r^(xfe-x) + rfp(x)iKdi 

+ d^mV^'^eix, - x) + rsP{x)\Vus\'' + snej|K.|'^(xe - x) 



+ d'mvcbmxt 



X] 



c^(x)[|v;,|^ + |v^. 



(16) 



For antineutrinos the virtual W is absorbed by positively charged quarks and antiquarks. 
In Eqs. ^ and ^ quantities like Vud are elements of the Cabibbo-Kobayashi-Maskawa 



(CKM) quark mixing matrix |2^. We have also introduced the so-called "slow rescaling'' 
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formalism P5|,PB| to account for threshold corrections in heavy quark production. For 
production of a heavy quark with current quark mass rrik, we define the quantities 




UQ') = 

Al ^ {MTn'-m% (17) 

In Eq. [l^, the quantity Mjp**^ is the minimum mass of the final state for the light quark 
to heavy quark transition. For the various quark flavors, we have Mjp*" ^ 2.8 GeV, 3 
GeV, and 6.2 GeV for transitions d — >• c, s — > c, and u ^ b, respectively (in this review 
we neglect all contributions from top quarks). 

Note that with this rescaling model, the Callan-Gross relation fails to apply in the 
region of heavy quark thresholds, even in the event that the interactions are purely 
transverse. The structure function ^ can be obtained from ^, Eq. |15|, by the 
replacement 

qi{x) -> 2xqi{x) q{^k) '^ikqi^k) , (18) 

and identical replacements for the antiquark distributions. With the same replacement, 
Eq. |TB|, one can obtain the structure function F^ ^ from F^ ^ . Similarly, the structure 
function \fY ^ can be obtained from F^ ^ by the replacement 

qk{x) -qk{x) (19) 



2.3 High energy limiting form for weak structure functions 

At sufficiently high momentum transfers, e.g. well above charm threshold, we have x. 
In the limit of very high momentum transfer, — *■ oo, we have ^fc(Q^) — *■ x and Xk —>■ 1 
for all heavy quark flavors. In this limit we once again recover the Callan-Gross relation 
(if (Jl/ctt 0). At very high momentum transfers and at very high energy, if we neglect 
correction terms of magnitude = |VJdP ~ sin(^c')^ 1 x 10""^, then the structure 

functions F^ (x, Q"^) reduce to 

lim F^''''{x,Q^) ^ dP{x)+uPix) + sP{x) + cP{x) 

Ff"^(x,g2) ^ uPix) + dP{x) + sP{x) + d'{x) 
^F^^P{x,Q^) ^ dP{x)-uP{x) + s''{x)-(^{x) 

^F^~''{x,Q'') ^ uP{x)-dP{x)-sP{x) + d'{x), (20) 

with the corresponding structure functions for neutrons obtained by replacing super- 
scripts p ^ n everywhere in Eqs. GUI Because of their simplicity we will generally use 
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Eq. ^ in deriving relations between structure functions, although we should revert to 
Eqs. |15| and |16| when comparing with data. This is particularly relevant for experiments 
at relatively low Q^, where threshold effects can be rather important. 
The assumption of charge symmetry for parton distributions is that 

= dP{x) 
^"■(x) = = s{x) 

c"(x) = c^ix) = c{x). (21) 

We have identical relations for antiquark distributions. With this assumption, all neutron 
parton distributions can be replaced by the corresponding distributions in the proton. To 
retain the charge symmetry violating parton distributions, we introduce the CSV parton 
distributions for up and down quarks via 

d'^ix) = uP{x) - 6u{x) 

m"(x) = d%x)-6d{x). (22) 

If the quantities Su{x) and 6d{x) vanish, then charge symmetry is exact. We have anal- 
ogous relations for CSV in antiquark distributions. We assume that the strange quark 
(and antiquark) distributions are the same in both the proton and neutron, as is given 
in Eq. ^ We make the same assumption for charm quarks. There is no theoretical or 
experimental reason to expect strange and charm distributions to vary from proton to 
neutron. 

It is useful to divide parton distributions into valence quark and sea quark parts. The 
valence up quark distribution in the proton is defined by u^{x) = u^{x) — u^{x). The 
valence quark distributions obey the following quark normalization conditions 

1 rl 

dxu^{x) = dx {u^{x) — u^{x)) = 
Jo 

^ dxd'^ix) = £ dx (rf"(x) - J"(x)) = 2; 

^ dxdP(x) = C (d^ix) - dP(x)) = 
Jo ^ ' 

dxulix) = f dx (u"(x) - u"(x)) = 1; 
Jo 

1 _ rl _ 

dx {s{x) — s{x)) = dx (c(x) — c(x)) = 0. (23) 









The CSV quantities defined in Eq. |22| can have both valence and sea pieces. The valence 
quark charge symmetry violating distributions are defined as 

6uv{x) = 6u{x) — 6u{x) 

6d,{x) = 6d{x)-6d{x) (24) 
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From these definitions of valence quark CSV, it is straiglitforward to sliow tliat tlie first 
moment of tlie valence quark CSV distributions (i.e., the integral over x) must vanish. 
We see that 



In Eq. ^ the integral over the valence quark distributions is fixed by the normalization 
condition to the number of valence up quarks in the proton (down quarks in the neutron). 
Since both of these are equal to 2, the integral of 6uv{x) must give zero. 

From Eq. ^ we see that the first moment of the heavy quark and antiquark distri- 
butions are identical. Until recently, it was customary to assume that the strange and 
charmed quark and antiquark distributions were equal for all values of x. That is, one 
assumed that 



with an identical relation for the charmed quark distributions. However, recently there 
has been both theoretical and experimental interest in whether the strange and antis- 
trange distributions are in fact equal. We will review how strange quark distributions 
are extracted, and the experimental situation regarding strange and antistrange quark 
distributions, in Sect. 2.6. 

2.4 CSV Contributions to Structure Functions 

In the high energy limit, well above heavy quark thresholds, the structure functions for 
charged current weak interactions on neutrons take the form 




(25) 



Six] = Six] = Six 



(26) 



lim Ff+"(x,g2) = d"(x) +m"(x) + s(x) +c(x) 



Ff~"(x,Q2) = m"(x) + c/"(x) + s(x) + c(x) 
^F3^^"(x,Q2) = d''{x)-u''{x) + s{x)-c{x) 

^F3^""(x,g2) = m"(x)- J"(x)-s(x) + c(x). 



(27) 



Introducing the CSV parton distributions from Eq. |2^, Eq. |2^ becomes 



lim Ff^"(x,g2) 

)2 — 5.^^ 



OO 



u^{x) + (F{x) + s(x) + c(x) 



5u{x) — 6d{x) 



dP{x) + uP{x) + s{x) + c{x) 
u^{x) — (F{x) + s{x) — c(x) 



dP{x) 



w 



(x) — s{x) + c(x) 



5d{x) — 5u{x) 
6u{x) + 6d{x) 



6d{x) + 6u{x) 



(28) 
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For completeness, we include the electromagnetic structure function for neutron tar- 
gets. Eq. in] becomes 



1/4 1 r - 

- - K(x) + + c{x) + c{x)] + - d"{x) + + s{x) + s{x) 

2 \t) 9 '- 



1 /4 



2 V9 
4 



dP{x) + dP{x) + c{x) + c{x)] + - y{x) + + s{x) + s{x)] 

9 



9 



6d{x) + 6d{x) — - [6u{x) + 6u{x)] 
9 



(29) 



Many tests of charge symmetry will involve deep inelastic scattering on isoscalar 
targets, which we label as iVg. Such reactions involve equal contributions from protons and 
neutrons. Under the assumptions we have listed previously, the weak and electromagnetic 
structure functions on isoscalar targets can be written (in terms of structure functions 
per nucleon) 



Fr^o(x,Q2) 
Ff-^o(a;,Q2) 
lFr''^ix,Q') 
1f3^-^o(x,Q^) 



^ (uP{x) + dP{x) + uPix) + dP{x) + 2s{x) + 2c{x) - 6u{x) - 6d{x)) 
i {uP{x) + dP{x) + uP{x) + (F(x) + 2s{x) + 2c(x) - 6d{x) - 6u{x)) 
i {uP{x) + d^{x) - vf{x) - d^ix) + 2s{x) - 2c{x) - 6u{x) + 6d{x)) 
i {uP{x) + dP{x) - u^'ix) - dF{x) - 2s{x) + 2c(x) - 5d{x) + 5u{x 



2 

1 /5 r 



4 V9 



u^{x) + d^{x) + u^{x) + d^{x) + - [s{x) + s{x)] 

9 



8 4 r - 1 1 

9 9 '- -'9 



(30) 



2.5 Isolating CSV Effects in Structure Functions 

In order to isolate and measure CSV effects, we need to find relations between various 
structure functions which depend on the validity of charge symmetry, and which can be 
tested. There are two such relations: the relation between the Fi (or F2) charge-changing 
electroweak structure functions from neutrino and antineutrino reactions, and the rela- 
tion between the F2 structure functions obtained from neutrino deep inelastic scattering, 
and the F2 structure function from deep inelastic scattering induced by charged leptons 
(muons or electrons). 

We first discuss the Fi structure functions for charge- changing weak interactions. For 
deep inelastic scattering on an isoscalar target, we can derive the following identity from 
Eq. H: 

hm (Ff^^o(x,g2)_irf-^o(a:,g2)) = hM,{^j:)-5u,{x)) 

+ s{x) — s{x) — c{x) + c{x) 

s[x) - s[x) - c[x) + c[x) (31) 
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For Eq. ^ to be valid, we must be at sufficiently high values of that we are well above 
both charm and bottom thresholds. Furthermore, we neglect terms of order = 
|Vtdp ^ 1 X 10^^. Eq. ^ should be true at all values of x. The final line of this equation 
holds in the limit that charge symmetry is exact. To avoid confusion in this review, we 
have introduced the notation *lf . This means that an equation is true provided charge 
symmetry is exact. In this way we hope one can distinguish between relations which 
are generally true, and those which require the (generally implicit) assumption of charge 
symmetry. 

At sufficiently high energies (well above heavy quark production thresholds) threshold 



effects should become negligible, and then Eq. |31| should be valid. If these structure 
functions are not equal at all values of x, this implies either charge symmetry violation in 
parton distributions, or inequality of the strange quark and antiquark distributions (the 
charm quark contributions should be quite small, and we know of no theoretical reason 
why charm and anticharm distributions should be unequal). In Sect. 5.1 we present 
theoretical estimates of valence quark CSV and strange/antistrange quark contributions 
to this relation. 

From Eqs. |2D| and ^ we can also derive relations between the F\ structure functions 
for neutrinos on protons, and antineutrinos on neutrons, 

''P(x, Q2) _ ""(x, Q2) = + ^u[x) + s[x) - six) - c[x) + c(x) 

^"(x, Q2) _ Ff''*'(x, Q2) = -bu[x) - 6d{x) + s{x) - s{x) - c{x) + c{x). (32) 

Eqs. ^ are valid under the same conditions as Eq. |31|, namely that we are well above 
heavy quark thresholds, and that we neglect Kobayashi-Maskawa matrix elements of 
order 10~^. If parton charge symmetry were exact, and strange quark and antiquark 
parton distributions are identical at all x, then we would expect 

F^'%x,Q^) = F^'-^ix^Q') (33) 

At sufficiently high energies, if charge symmetry is valid for valence quark distribu- 
tions, and if the strange quark and antiquark distributions are equal at all x, then the 

Fr 

structure functions are identical for isoscalar nuclear targets. Identical relations 
hold for the F2 structure functions, when we include the longitudinal/transverse ratio 
R of Eq. |T^. These equations have been used to extract the structure functions F3 in 
electroweak reactions. Using Eqs. ^ and ^, we can derive 

3n 



2G^MmE 



{da^P/dx - da^'^'/dx) = F^ ^(x, Q') - Ff ""(x, Q 

1 



2G^MnE 



-r -y.Fr'ix,Q')+xFr''{x,Q' 

I' l{xFr'{x,Q')+xFr''{x,Q'- 
{da'^'/dx - da^P/dx) = Ff ^"(x, Q^) - Ff "^(x, Q^) 
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+ 

cs 



\{xFr''{x,Q')+xFr'{x,Q'] 

UxFr''{x,Q')+xFr'{x,Q' 



3tt 



2G'^MnE 



[da''''' /dx - da'""" /dx) = Ff"-^" (a;, Q^) _ f]^-a^o g2 



cs 



where in Eq. ^ we define 

da''P{x,Q^)/dx 



dy d'^a'^^^x, y)/dxdy 



,(34) 



2M^ sin^ e 



(35) 



From Eqs. and we see that if charge symmetry is exact, and if the strange quark 
and antiquark distributions are equal at all x, then by taking the difference between cross 
sections for the appropriate charged current cross sections for neutrinos and antineutrinos, 
the relevant F2 structure functions will cancel, leaving just the F3 structure functions. 
This follows from Eq. In particular, the difference between the charged current 
cross section from neutrino scattering on an isoscalar target, and the cross section from 
antineutrinos on that target, is just equal to the sum of the structure functions xF^ for 
neutrinos and antineutrinos on the isoscalar target. 

If we don't require charge symmetry and equality of strange and antistrange quark 
distributions, then Eq. becomes 



2G^MnE 



3tt 



2G^MnE 



{da'^P/dx - da^''/dx) 



(da'^^'/dx-da^P/dx) 



2 

2x 
1 



2G^MnE 



[da^^'/dx-da^'^ydx 



2x [6d{x) + 6u{x) + s{x) — s{x)] 

l{xFr'ix,Q') + xFr''{x,Q' 
—6u{x) — 6d{x) + s{x) — s{x) 

'xFr'\x,Q')+xFr'{x,Q'^ 
= X [6d^{x) — 5u^{x) + 2 (s(x) — s(x))] 
+ ^(xF3^'^^"(x,Q^)+xF3"(a:,Q^)), (36) 



There are additional terms arising from charm quark distributions; they are given in Eq. 
32| , but we have not included them in Eq. 

From Eq. ^ we see that the difference between neutrino and antineutrino cross sec- 
tions on an isoscalar target contains not only the F^ structure functions, but two residual 
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terms - one of which depends on the quark CSV amplitudes, and the other depending 
on the difference between strange and antistrange quark distributions. These additional 
terms then make a contribution at various x values to the sum of the F3 structure func- 
tions on isoscalar targets. From Eq. without these additional terms the difference 
between charge-changing reactions induced by neutrinos and antineutrinos gives the sum 
of up plus down valence quark distributions in the nucleon. 

The other relation between structure functions which allows an experimental test of 
charge symmetry is the so-called "charge ratio" or the "5/18*'^ rule." Neglecting for the 
moment the longitudinal/transverse ratio R (which will cancel if we form the ratio of the 
two structure functions), we have from Eq. ^ 



F2 [x, Q 



X 



Q{x) 



Q{x) ^ 



- (6u{x) + 6u{x) + 6d{x) + 6d{x] 
1 



X 



6 



{s{x) + s{x) — c{x) — c{x)) 



2 1 

- (5d{x) + 6d{x)) — — {6u{x) + 5u{x)) 
9 ^ ^ 18 

gj{x) + qj{x) 

j=u,d,s,c 



(37) 



From Eq. p^, if we take the ratio of the two F2 structure functions we obtain 



1 - 



3 {s{x) + s{x)) A5d{x) + A5d{x) + 5u{x) + 5u{x) 
5Q(s) hQ{x) 

5d{x) + 5d{x) + 5u{x) + 5u{x) 



2Qix) 

^ 3 (s(x) + s{x)) ^ 3 {6u{x) + 6u{x) - 6d{x) - 6d{x) 



cs 
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1 - 



5Qix) 

3 (s(x) -|- s{x)) 
5Q{x) 



WQ{x) 



{31 



In Eq. ^ we have dropped the small contribution from charm quarks. In the second 
of Eqs. ^ we have expanded to lowest order in the (small) CSV terms. The ratio 5/18 
in Eq. ^ the so-called "charge ratio" for these structure functions, occurs because the 
virtual photon couples to the squared charge of the quarks, while the charged-current 
reactions induced by neutrinos couple to the weak isospin mediated by W exchange. 

In the charge symmetric limit, we can use the F2 structure functions from either 
neutrino or antineutrino induced reactions, since 



17 



If we use the neutrino or antineutrino F2 structure functions instead of their average in 
Eq. the only thing which changes in the ratio is the weighting of the various CSV 
terms, plus an additional contribution if the strange and antistrange distributions are not 
identical. 

In order to get theoretical estimates for the structure function relations, we need 
to know the quark CSV contributions, and we also need to know the strange quark and 
antiquark distributions in the nucleon. The strange quark and antiquark distributions can 
be obtained experimentally by measuring production of opposite sign dimuons in reactions 
induced by neutrinos or antineutrinos. We review this in the following section. We 
obtain estimates for the charge-symmetry violating parton distributions using the model 
for quark CSV which we derive in Sect. 3. In Sect. 4, we will estimate the contribution 
of the CSV terms to existing tests of charge symmetry at high energies. In Sect. 6, we 
will discuss possible CSV sea quark contributions to the various sum rules, in particular 
the Adler and Gross-Llewellyn Smith sum rules. 



2.6 Extraction of Strange Quark Distributions 

In Sect. 2.3 we introduced the assumption that strange quark and antiquark distributions 
were identical at all x, and were identical for protons and neutrons (a similar assumption 
is made for charmed quark distributions). In Sect. 2.5 we saw that structure function 
tests of charge symmetry also contained contributions from strange quark and antiquark 
distributions. Thus, to get accurate tests of parton charge symmetry, we must have 
reliable measurements of strange quark distributions. 

The strange quark distribution can be assessed in two ways. First, it can be obtained 
"indirectly" by comparing DIS reactions induced by charged leptons (muons or electrons) 
with charge- changing currents from neutrinos. From Eq. |3^, we see that for an isoscalar 
target (assuming the Callan-Gross relation and neglecting charm quark contributions), 
we have 

6u{x) + 6u{x) — 6d{x) — 6d{x) 
+ ^[six)+s{x)] (39) 

Eq. ^ is "indirect" , in that we must compare experiments with muon and neutrino beams, 
performed under different conditions and with different normalizations. Furthermore, we 
need to know the CSV contributions in order to extract the strange quark distribution. 

The "direct" way of extracting strange quark distributions is by measuring the yield 
of opposite sign dimuons produced in nuclear reactions induced by neutrinos. In leading 
order, the incoming neutrino has a hard scattering with an s or ci quark, producing a 
charm quark which fragments into a charmed hadron. The subsequent semileptonic decay 
of the charmed hadron produces an opposite sign muon, through the process 

Uf, + N ^ fi' + c + X 



F-,''''\x,Q')-^Fr\x,Q') = ^ 
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^ )= )x 



Figure 4: Contribution to production of an opposite sign dimuon pair via intermediate 
charm quark production. Leading-order quark initiated diagram, plus two NLO gluon- 
initiated diagrams. 



The antineutrino process produces charm antiquarks from d and s sea quarks. As the first 
muon produced tends to come from the original scattering and to have larger transverse 
momentum with respect to the direction of the hadron shower, experimenters can tell 
whether the muon pair arose from a neutrino or antineutrino collision. In this section we 
summarize the latest direct measurements and their results. 

In Fig. ^ we show the leading order [LO] mechanism for production of ^~ pairs in 
neutrino-induced reactions. A virtual is absorbed on an s or quark, producing a 
charm quark, which then decays semi-leptonically. Using the slow rescaling model 
and assuming charge symmetry, the leading order cross section for production of opposite 
sign dimuons, for neutrino reactions on an isoscalar target, has the form 



2 

mi 



D{z)B, (40) 



In Eq. D{z) is the fragmentation function for a charmed quark into a charmed hadron, 
and Be is the branching ratio for semileptonic decay of a charmed hadron. The quantity 
is given by the slow rescaling formula, Eq. [l^, and the parton distributions depend 
through QCD on the scale /x^. The corresponding cross sections for antineutrino inter- 
actions are obtained by replacing qki^c^ 1^) by qkiic^ 1^) for each quark flavor, in Eq. 

m 

Because the CKM matrix element Vcs is substantially larger than Vcd-, i-e. \Vcs\^ ~ 0.95 
while \Vcd\^ ~ 0.05, opposite sign dimuon production from neutrinos is most sensitive to 
the strange quark distribution in the nucleon, even though the d quark content of the 
proton is roughly ten times the strange quark density. So the Cabibbo suppression of 
the d quark contribution to charm production makes the strange quark contribution 
relatively more important. This suppression factor is also present for reactions induced 



19 





Figure 5: a) Sea quark distribution xJ2j 3 = u, d, s, from NLO analysis of CCFR 

Collaboration data, Ref. [BO], b) Strange quark distribution s{x) from same data. 



by antineutrinos, but here the relevant antiquark distributions s and d are equal to within 
about a factor of two. 

Because of the suppression of contributions from the valence quarks, the next-to- 



leading [NLO] contributions from gluon exchange turn out to be quite important |27 



The most important such processes, t and u channel diagrams initiated by gluons, are 
also shown in Fig. ^ ^ . The extra factor of as which enters these diagrams is compensated 
for by the fact that the gluon density is an order of magnitude larger than the antiquark 
density. It turns out to be crucial to include the NLO contributions to this reaction, in 
this kinematic region (relatively near charm threshold). 

Recent measurements by the CCFR collaboration (experiments E744 and E770 at the 
Fermilab Tevatron with the quadrupole triplet neutrino beam) [^T],^^ have amassed a 
substantial number of events for both neutrinos and antineutrinos. The following results 
have been obtained with the latest NLO analysis of the CCFR data In Fig. §a we 
plot xq{x, /i^) = 
for a scale /x^ - 
the same scale 



xu{x, fi ) + xd{x, fi ) +xs{x, n ) extracted from these measurements 



at 



4 GeV^. In Fig. |]b we plot the strange quark distribution xs{x,fi 
Both the LO results, and the NLO results, are plotted. We see that the 
NLO curves differ dramatically from the LO results, with the NLO curves much softer 
than the LO results (they are significantly larger at small x and fall off faster) Q 
Second, the strange/nonstrange antiquark fraction can be extracted. If we define 



S + S 



where 



^We have not shown radiative-gluon and self-energy diagrams which also occur in NLO. 
■^However, see the paper by Gliick et al. ||3l[| , who claim that a consistent treatment of acceptance 
corrections gives NLO results, for the same CCFR data, which are much closer to the original LO results. 
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S = xs(x)dx, the CCFR group obtain 

JO 

,„„+0.046, , +0.023, , , , 

K = 0.477 _ Stat + fsyst 41 

-0.044^ ^ -0.024^ ^ ^ ^ ^ 

The extracted value of k shows a substantial violation of SU(3) symmetry in the 
nucleon sea. However, the value of k obtained in the NLO analysis |]2T| is substantially 



larger than the result obtained in a LO analysis of the same data I^Sf, which produced a 
value ~ 0.37. This arises because the CCFR structure functions extracted from the 
NLO analysis give larger values at small x for both the nonstrange sea and the strange 
sea. The importance of the NLO analysis is quite striking here. 

Next, the CCFR group compared the x dependence of the strange and nonstrange 
sea. They parameterized the strange quark distribution as 



xs{x, /i^) = ^^(1 — x] 



a 



xu{x, /i^) + xd{x, /i^ 



(42) 



They obtained the best fit value a = —0.02 + ^^'^^(stat) + ^.oll i^Y^^)- '^^^^ value is 
consistent with zero, so the x dependence of strange and nonstrange sea appears to be 
identical. The LO analysis |^ appeared to show a much softer s distribution than the 



nonstrange sea. 

Fourth, this data can be used to extract the CKM matrix element Vcd- At present, 
the uncertainty in this parameter gives the greatest contribution to the uncertainty in 
the Weinberg angle, sin^(6'vi/)- 

Finally, the CCFR group tested whether the strange antiquark distribution s{x) dif- 
fered in shape from the strange quark distribution s{x). The LO analysis of this group 
had suggested some difference between the two distributions. There have been several 
theoretical suggestions that this might be the case [p^-|36[|. If this should prove to be 



the case, then our formulae for CSV need to be modified to include differences between 
quark and antiquark distributions for heavy quark flavors. This will turn out to be par- 
ticularly important in tests of charge symmetry involving charge- changing weak currents 
on isoscalar nuclei, as is discussed in Sect. 5.1. 

The CCFR group analyzed their strange quark distributions assuming that xs{x, /i^) = 
(1 - x)'^xs{x,^i^). They obtained the value (3 = -0.46 ± 0.42 ± 0.36 ± 0.65 ± 0.17; the 
quoted errors (from left to right) are statistical, systematic, from the uncertainty in 
the semileptonic charged hadron branching ratios, and the uncertainty in /i^ scale. The 
value obtained is consistent with zero, i.e. no difference in the shape of the strange and 
antistrange distributions. 

Let us review the outline for subsequent sections of this paper. In Sect. 3, we review 
model calculations which give order of magnitude estimates for charge symmetry violating 
contributions to both valence and sea quark distributions. In Sect. 4 we review existing 
experiments which allow us to put upper limits on the magnitude of charge symmetry 
violation in the parton distributions. We will show that present upper limits on CSV in 
parton distributions depend on the x region in question. In the region 0.1 < a; < 0.4, 
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the comparison of CCFR neutrino structure functions to the structure functions from the 
NMC muon measurements suggests an upper hmit on CSV of a few percent. At larger 
values of x, there is much less experimental data, so the upper limits on parton CSV are 
at least 10% in the parton distributions. For smaller values 0.01 < x < 0.1, the NMC 
and CCFR results currently disagree at a level between 10 — 20%. 

We will also review recent experimental tests of flavor symmetry in the proton sea. 
We will show that these "tests" demonstrate large breaking of either flavor symmetry or 
charge symmetry, and that at the present time it is difficult to rule out significant breaking 
of charge symmetry in the nucleon sea (although we certainly do not expect CSV effects 
of the magnitude necessary to agree with recent Drell-Yan experiments). In Sect. 5 
we suggest several new experiments which can differentiate between charge symmetry 
and flavor symmetry. These experiments, if carried out, could either demonstrate the 
existence of CSV terms in parton distributions, or provide more stringent upper limits 
on quark CSV. 

3 Theoretical Estimates of Charge Symmetry Viola- 
tion in Parton Distributions 



In this section, we will derive theoretical estimates for charge symmetry violation in 
parton distributions. First, we will review theoretical calculations of CSV for valence 
quark distributions. Next, we will give an estimate of the size of CSV effects in antiquark 
distributions. We will discuss the robustness of these estimates. In later sections, these 
calculations will be used to provide estimates of the magnitude of CSV effects which 
could be expected in various experiments. In absence of more detailed calculations of 
CSV effects, and lacking flrm upper limits from experiment, these estimates are the best 
we have at present. As we will show, we believe that our estimates of both the sign and 
magnitude of valence quark CSV effects should be rather well determined. Neither the 
sign nor magnitude of sea quark CSV contributions is well known; however, our model 
calculations predict very small CSV effects from the sea. 



3.1 Charge Symmetry Violation for Valence Quarks 
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Theoretical investigations of parton CSV for valence quark distributions by Sather 
Rodionov, Thomas and Londergan |T^ and Benesh and Goldman concluded that 
one could make reasonably model-independent estimates of the size of these effects. Here 
we follow the method for calculating twist-two parton distributions with proper support, 

MM 



which has been developed by the Adelaide group 
evaluation of quark distributions through the relation 



The starting point is the 



q{x,^i') = MJ2 I < ^l^+(0)|iV > \'6{M{l-x)-p^) 



(43) 



X 
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In Eq. ^ ip^ = {1 + a3)4'/'2, X represents a complete set of eigenstates of the Hamilto- 
nian, and /z^ represents the starting scale for the quark distribution. 

The advantage of this method is that the resulting parton distribution /x^) is 
guaranteed to have proper support, i.e. it vanishes for x > 1, regardless of the model 
used for the matrix element < X\ip^{0)\N > in Eq. 43. Parton distributions calculated 
from quark models generally lack this support, and this can lead to serious problems 
in obtaining reliable results. Thus, Thomas and collaborators showed that reasonable 



parton distributions could be obtained from models such as the MIT bag |BB[ . The other 



advantage of this method is that it is often possible to obtain reasonable results taking 
into account only the lowest-energy spectator states in the sum over states of Eq. |43. 

We want to use the relation in Eq. ^ to calculate differences in parton distributions 
due to violation of charge symmetry. Thus, we wish to estimate the difference between, 
say, the up quark distribution in the proton and the down quark distribution in the 
neutron. From Eq. ^ we see that CSV contributions will have four sources: first, from 
electromagnetic effects which break charge symmetry; second, from charge symmetry vi- 
olating mass differences of the struck quark; third, from mass differences in the spectator 
multiquark system; and fourth from charge symmetry violation in the quark wavefunc- 
tions. In model calculations, it was found that the quark wavefunctions are almost 
invariant under small mass changes typical of CSV. At high energies, electromagnetic ef- 
fects should also be small, and we neglect these. Consequently, parton charge symmetry 
violation will arise predominantly through mass differences of the struck quark, and from 
mass differences in the spectator multiquark system. 

As an example, we consider valence quark CSV where for the intermediate states 
X in Eq. ^ we include the lowest two-quark spectator states from the MIT bag model 
||T9|, |38|, |39|, ^ . There are more sophisticated quark models available but the similarity 



of the results obtained by Naar and Birse using the color dielectric model suggests 
that similar results would be obtained in any relativistic model based on confined current 
quarks. In Fig. ^we show schematically the lowest contributions to the "majority valence 
quark" distributions, i.e. and d'^{x). The majority quark CSV term is as defined 

in Eq. 6uv{x) = u''^{x) — From Fig. ^ we see that the only contribution to 

the "majority" quark CSV is the up-down mass difference m„ — rrid', the intermediate 
spectator diquark is the same [ud) for both proton and neutron. 

In Fig. 1^ we show the lowest contributions to the "minority valence quark" distri- 
butions, i.e. d^{x) and From Fig. ^ we see that there are two contributions to 
minority quark CSV. One is the up-down mass difference m„ — m^; the second source of 
charge symmetry violation comes from the intermediate spectator diquark mass, which 
is uu for the proton and dd for the neutron. In Fig. |] we show the calculated minority 
valence quark CSV term, 6d^{x) = d^{x) — u'^{x). Fig. ^ shows the majority valence 
quark CSV term, 6uv{x) = u^{x) — d^i^x). These are calculated from Eq. ^ using quark 
wave functions from the MIT bag model. The contributions are calculated at the bag 
scale, then evolved upwards in Q^. At small x, 6d^{x) is negative, while for larger x it 
is positive. The majority quark CSV term has exactly the opposite sign. As a result. 



23 



Figure 6: Schematic picture of lowest energy contributions to "majority" quark distribu- 
tions: struck quark plus spectator diquark. (a) for proton, (b) d"'{x) for neutron. 




Figure 7: Schematic picture of lowest energy contributions to "minority" quark distribu- 
tions, (a) d'^{x) for proton, (b) u"'{x) for neutron. 
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Figure 8: Calculation of minority valence quark CSV contribution, Sdjj{x) vs. x. Dotted 
curve: = iq QeV^; dashed curve: = 100 GeV^; dot-dashed curve: = 400 GeV^; 
dash-triple dot curve: = 10,000 GeV^. From Ref. [|7|. 



theoretical calculations suggest that 



\6uv{x) + 6d^{x)\ « \Suv{x) — Sdv{x)\ 



(44) 



so that experimental quantities which depend on the sum of the majority and minority 
valence quark CSV terms should be substantially smaller than those which depend on 
the difference between the majority and minority CSV terms. 

Because the integral over x of the valence quark distributions is normalized (one for 
the minority valence quark distribution, two for the majority distribution), the integral 
of the CSV distributions must be zero, i.e. 



1 rl 

6uv(x) dx = 6dJx) dx = ^ 
Jo 



(45) 



In Fig. |1^ we show the percent CSV contribution, 6d^{x)/d^{x) vs. x. For large x, i.e. 
X > 0.5, the minority valence quark CSV is predicted to be between 5 — 10%. This is 
an extremely large violation of charge symmetry, especially since at low energy scales 
(e.g., low energy nuclear physics) charge symmetry is generally valid to at least 1%. 
Compared with the minority quark CSV contributions, the majority quark term has 
roughly the same magnitude, as can be seen from Figs. § and |^. However, for large x the 
percent majority quark CSV is predicted to be much smaller than the minority quark 
CSV fraction. 

Why should the minority valence quark CSV, 6dv{x), be a much larger percentage 
at large x than the majority CSV term, Suv{x)7 We can understand this qualitatively 
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Figure 9: Majority valence quark CSV contribution, 6uy{x). Notation is that of Fig. ||. 
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Figure 10: Model calculation of percent minority valence quark CSV term, 6dy{x)/dy{x) 
vs. X, evolved to = 10 GeV^. From Ref. 0. 
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because the dominant source of CSV is the mass difference of the residual diquark pair 
when one quark is hit in the deep-inelastic process. For the majority quark distribu- 
tion, the diquark pair is ud for both proton and neutron, as can be seen from Fig. ^. 
Consequently there is no CSV contribution from this term. For the minority quark dis- 
tribution the residual diquark is uu in the proton, and dd in the neutron. Thus, in the 
difference, d^{x)—u'^{x), the up-down mass difference enters twice. Including electromag- 
netic effects, we know that the diquark mass difference is SM2 = Mf - M^"" ^ 4 MeV. 
From Eq. ^ we can see that the valence quark distribution q^lx) will peak at roughly 
Xpeak ~ 1 — M2/M where M is the nucleon mass. For the minority quark distribution, 
say, (iv(a;), the lowest diquark state for the uu pair will be = 1 by the Pauli principle. 
From the — A mass splitting we know that an = 1 diquark pair will have an effective 
mass M2{S = 1) ~ 800 MeV, while the S = diquark (which is present only for the 
majority quark term u^{x)) has a mass M2{S = 0) ~ 600 MeV. 

As was pointed out by Close and Thomas |^3|, this explains why the down quark 
distribution d^{x) peaks at a value x ~ 0.1, while the up valence quark distribution 
u^{x) peaks at x ~ 0.3. It also predicts that d^{x) << u^{x) at large X, cLS IS observed 
experimentally. We can obtain an estimate of the magnitude of minority quark charge 
symmetry breaking, 

Sd^ ^ 6M2 

dl{x) ^ M-M2 ^ ' 

Eq. ^ gives an estimate of about 4% for minority quark CSV. This can be compared 
with the 3-7% estimates of CSV by Sather |jl8|, Rodionov et al. [jl9l, and Benesh and 



Goldman 



Fig. [T^ shows the fractional minority quark CSV term, 2{d^ — u'l)/{d^ + -u") vs. x 
for several values of the intermediate mean diquark mass. Although the precise value of 
the minority quark CSV changes with mean diquark mass, the size is always roughly the 
same and the sign is unchanged. This shows that "smearing" the mean diquark mass 
will not dramatically diminish the magnitude of the minority quark CSV term (the mean 
diquark mass must be roughly 800 MeV in the S = 1 state to give the correct A^ — A 
mass splitting). 

We reiterate that the magnitude of charge symmetry violation, predicted for the 
"minority" valence quarks at large x, is both extremely large and surprising. First, since 
the integral over x of the valence quark CSV term is zero, as given in Eq. large CSV 
contributions to valence quark parton distributions do not necessarily imply large charge 
symmetry violation at low energies. Second, experimental verification of these CSV 
effects is not a simple matter. Since we predict that charge symmetry is well obeyed for 
the majority valence quarks, finding CSV effects requires experiments which are sensitive 
to the minority quark distributions, in a region (large Bjorken x) where the minority 
quark distributions are much smaller than the majority quark terms. In Sects. 4 and 5, 
we will review the current experimental limits on parton charge symmetry, and suggest 
experiments which would be sensitive to our predicted effects. 

However, as we have tried to emphasize, our predictions of charge symmetry violation 
depend on rather simple assumptions, which have been shown to produce reasonable 
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parton distributions Furthermore, all theoretical calculations of parton charge 



symmetry obtain the same qualitative result: that the percentage of CSV for minority 
valence quarks should be substantially larger than the fraction of majority quark CSV 
at large x [^,|19|,|3^. To state this another way, theory predicts that 6d^{x) ~ —6u^{x); 
however, at large x d^{x) « Uy{x), so the fraction of minority quark CSV will be 
substantially larger than for majority quarks. In particular, the calculation of Sather [ P^ 
gave what was essentially a model-independent estimate of parton CSV. Sather also 
predicted that minority quark CSV effects should be, relatively, substantially larger than 
majority quark CSV. This gives us confidence that our predictions are robust. 

Finally, since all phenomenological fits to date assume parton charge symmetry at 
the outset, existing parton distributions then have CSV effects included implicitly. As 
we pointed out in Sect. 2, a truly consistent treatment of these effects would begin at the 
outset with charge asymmetric parton distributions and proceed to fit experimental data, 
without making the prior assumption of charge symmetry. Then the CSV contributions 
could be deduced in a consistent fashion. In the absence of such a consistent procedure, 
our calculations of charge symmetry violation should be taken as order of magnitude 
estimates only. This is a difficult process, and we assume it will not take place until 
there is some definite experimental evidence of charge symmetry violation in parton 
distributions. 



3.2 Estimate of Charge Symmetry Violation for Sea Quarks 

In the preceding section, we made estimates of valence quark charge symmetry violation, 
using a formalism for quark distributions which was guaranteed to produce the proper 
support. We could use the same formalism to calculate antiquark distributions, 

q{x, fi') = MY.\<X\^pi{0)\N>\' 5{M{1 + x) - p+ ) (47) 

X 

We would proceed with the same assumptions as for valence quarks: we would take light 
cone quark momentum wave functions from simple models, and truncate the sum over 



states in Eq. to the lowest energy states which can contribute. 



Recently, such a procedure has been followed by Benesh and Londergan |Q. One is 
confident that reasonably model-independent estimates can be made for valence quark 
CSV. This is borne out by substantial agreement between various theoretical results. 
Calculations of sea quark CSV require some additional assumptions, and the model de- 
pendence of sea quark CSV is not entirely clear. However, in these calculations one 
predicts that charge symmetry violation for antiquarks should be very small, probably 
at least an order of magnitude smaller than the corresponding effects for valence quarks. 

We can give simple qualitative arguments why sea quark CSV effects should be small. 
The relative magnitude of CSV effects will be given approximately by 

q M 
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where M is the energy of the lowest intermediate state in Eq. and 5M is the mass 
difference for intermediate states related by charge symmetry. For antiquarks, the lowest 
energy states which contribute are four-quark states, whose energy is roughly twice the 
energy of the lowest diquark states which contribute for the valence quarks. We estimate 
the mass difference between charge symmetric four quark states as 5M ^1.3 MeV, or 
three times smaller than the mass difference for minority valence quarks. Thus, a naive 
expectation would be that CSV effects for sea quarks would be at least a factor of six 
smaller than for minority valence quarks. 

In their calculations Benesh and Londergan typically obtained estimates for 
sea quark CSV at least an order of magnitude smaller than those for valence quark 
CSV effects. As a result, in subsequent sections we will occasionally neglect sea quark 
CSV effects relative to valence quark CSV terms. Benesh and Londergan also obtain 
5u{x) ~ — (5(J(x); therefore, one would expect that observables proportional to the dif- 
ference between up and down sea quark CSV distributions would be significantly larger 
than observables which measure the sum of sea quark CSV terms. In Sects. 6 and 7 we 
review parton sum rules. There we will show that, in principle, it would be possible to get 
explicit measurements of charge symmetry and/or flavor symmetry violation in parton 
distributions, by measuring appropriate integrals of deep inelastic cross sections. 



4 Present Experimental Limits on Parton Flavor Sym- 
metry and Charge Symmetry 

In this section, we review the limits we can place on parton flavor symmetry in the 
proton sea, and on charge symmetry of parton distributions. In the following section, 
we propose a series of experiments which could in principle sharpen the limits on CSV 
in parton distribution, and which could also discriminate between flavor symmetry and 
charge symmetry violation. 

4.1 Drell-Yan Tests of Flavor Symmetry in the Proton Sea 

Recently, there has been much interest in the details of proton sea quark distributions. 
The NMC experiment measured the F2 structure functions for muons on proton and 

deuteron targets. With these data they were able to test the Gottfried Sum Rule , which 
requires the integral of the difference between proton and neutron F2 structure functions. 
The experimental value Sq = 0.235 ± 0.026 is more than four standard deviations below 
the "naive" expectation of 1/3. Both the experimental and theoretical situations are 
summarized in detail in the recent review by Kumano [^]. We review the Gottfried Sum 
Rule in Sect. 6.1. 

The most likely cause for the NMC result is a substantial difference in the d and u 
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distributions in the proton. The NMC experiment suggests that 







dx 



0.147 ±0.039 



There is no leading order QCD correction to the Gottfried Sum Rule. Ross and Sachrajda 



47| showed that higher order perturbative QCD contributions lead to a value much 
smaller than this. Consequently, there has been much interest in experiments which 
might give a "direct" measurement of the sea quark distributions dP{x) and vP{x), and 
which might map out their x dependence (the NMC experiment gives only the integral 
over X of this difference). 

Ellis and Stirling |]l3l suggested that this could be measured by comparing Drell-Yan 
processes initiated by protons, on proton and deuteron targets. We review here the 
information which could be obtained from these measurements. 

In the Drell-Yan [DY] process |^ one observes lepton pairs of opposite charge and 
large invariant mass which arise from hadronic collisions. This process occurs when a 
quark (antiquark) from the projectile annihilates an antiquark (quark) of the same flavor 
from the target. This produces a virtual photon which subsequently decays into a pair of 
charged leptons. The process is shown schematically in Fig. [TTla, for NN DY processes. 
A quark in one nucleon annihilates an antiquark of the same flavor in the other nucleon. 
In Fig. p!T| b we show the corresponding DY process for vr^ + p reactions, in the valence 
region for both particles. 

The Drell-Yan process for the interaction of hadron A with hadron B has the form 



rfV^^ Anal 



dxi dx2 9sxi X2 



K{x,, X2) V ei qt{xr)q^{x2) + qt{x,)qf{x,) (49) 



In Eq. s is the square of the CM energy, and xi and X2 are, respectively, the longitudi- 
nal momentum fractions carried by the target (projectile) quarks (or antiquarks) of flavor 
i and charge Cj. For example, the quantity qf{x2) is the antiquark distribution of the 
target for quarks of flavor i and momentum fraction X2- The factor K{xi,X2) accounts 
for the higher-order QCD corrections which enter the DY process. Detailed reviews of 



their form can be found in several articles [pO|,pT 



The values of Xi and X2 can be extracted from experiment through the equations 

= sxiX2^2Pi+Pi^{l-cos9i,+(,-) 
2(P,+ +P,-)^ 

Xp = Xi — X2 = — 1 

S 

T = X1X2 

In {X1/X2) , . 

y = — 2 — ^ ' 

In Eq. Pi+ and Pg- are respectively the laboratory momenta of the outgoing leptons, 
{Pi+ +Pi-)j^ is the longitudinal momentum of the lepton pair, and 9i+i- is the angle 
between their momentum vectors. 
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Figure 11: Schematic picture of Drell-Yan [DY] process for two hadrons. (a) A^A^ DY is 
sensitive to antiquarks in one nucleon; (b) — p DY process, in the valence region for 
both particles. 



The experiments have compared Drell-Yan cross sections for incident protons on pro- 
ton and deuterium targets. Taking ratios of Drell-Yan cross sections avoids the necessity 
for precise knowledge of the factor K in Eq. Assuming the validity of the impulse 
approximation, the DY cross section on the deuteron is just the sum of the DY cross sec- 
tions on the free proton and neutron. In that case the DY cross sections are proportional 
to 



a 



pp 

DY 



pD 
<^DY 



^ [U^{XI)U'{X2) + tP{x{)u\x2)] + \ \(F{xi)d''{x2) + $'{xi)(F{x2) 

9 9 

^ [uP{xi) {UP{X2) + m"(x2)) + K(x2) + ^"(xa))] 

9 

1 r 
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<ff\xi) (cF(x2) + J"(a;2) ) J^(xi) {<F{x2) + ^"(xs)) 



(51) 



If we assume charge symmetry then Eq. ^ reduces to 



a 



pD 
DY 



cs 



+ 



uP(x,) + l^dP{x,)) {uP{x2)+dPix2)) 
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1 ^„ 



^U^ix,) + -dPix,)] iuP{x2) + dPix2)) 



(52) 



The physics is most clear if we go to large Xi, i.e. large x for the projectile proton, 
and substantially smaller X2 for the target. In this regime the probability for finding 
antiquarks in the projectile is extremely small, so as a good approximation the DY 
process proceeds by quarks in the projectile annihilating antiquarks in the target. In this 
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region the ratio of DY processes on the deuteron, to those from the proton, are given by 
B (. cs + (l+r(.0/4) 



In Eq. ^ we define the ratios 

= ^ 

r(^) = ^ (54) 

From Eqs. ^ and ^ we see that if dP{x) = vF{x) for some x = Xe, then Rsea{xi,Xe) = 
1. For large x, the ratio r(a;) is small (the probability of finding an up quark in the proton 
at large x is significantly higher than for finding a down quark). In that case, as Xi —>■ 1, 
we have 

Rsea{Xl,X2) Z ^ {l + R{X2)) (^l + [l - R{x2)]^ , (55) 

From Eq. we see that the ratio of DY cross sections in this kinematic region should 
directly measure the ratio of the down antiquark to up antiquark distributions in the 
proton, at a given value of X2- Neglecting terms of order r{xi), the quantity Rsea{xi,X2) 
would be less than one if dP{x2) < vP{x2), and greater than one if dP{x2) > vP{x2). There 
have been three recent experiments which enable us to extract the difference between 
down and up antiquark distributions in the proton. 



The experiment E772 at FNAL measured Drell-Yan processes for 800 GeV protons 
on a variety of nuclear targets. The targets employed were D, C, Ca, Fe, and W. As there 
was no proton target, the D /p ratio could be inferred by comparing the isoscalar targets 
with those for which N ^ Z. For non-isoscalar targets the excess neutron fraction is 
proportional to e = {N — Z)/A. As e is rather small for these targets, it is difficult to 
measure neutron/proton differences, and hence hard to isolate any difference between 
down and up antiquarks in the proton. The experimental results appeared to show 



differences between dP and vP [^, and did disagree with some theoretical suggestions 
for dF /vF, but it was hard to draw any firm conclusions regarding this question from the 
E772 experiment. 



Experiment NA51 at CERN measured Drell-Yan processes for 450 GeV protons 
on proton and deuteron targets. The NA51 data looked primarily at symmetric kine- 
matics Xi = X2- The symmetric geometry is particularly good for minimizing errors in 
comparison of different experiments. For this geometry the approximations used to gen- 
erate Eq. ^ are not valid. Instead, Ellis and Stirling |13[ showed that for symmetric 



kinematics the ratio of Drell-Yan cross sections could be written as 



pD 



x,x) _ 1 / l{R{x)+r{x)] 
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Figure 12: The ratio of pD to pp DY cross sections as a function of target X2, preliminary 
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results of the E866 Collaboration 
distributions. 



Curves are for phenomenological parton 



The NA51 group obtained a ratio for a single averaged point < x >= 0.18. From 
their measured asymmetry in Drell-Yan processes, they extracted the value 



dP 



|<a;>=0.18 



0.51 ± OM{stat) ± 0.05{syst) 



(57) 



The NA51 result, although for only a single average x value, suggests that there are twice 
as many down antiquarks as up antiquarks at x = 0.18. 

The E866 group at FNAL |jl5| compared Drell-Yan processes for 800 GeV protons 
on liquid hydrogen and Deuterium targets. The E866 experiment has both the high 
statistics and the wide kinematic range which made it difficult for prior experiments to 
investigate this issue. In Fig. |T2| we show preliminary results from E866 for the ratio 
of DY cross sections a^^ /2a^^, for positive xp = xi — X2- In this kinematic region we 
expect the antiquarks to come predominantly from the target. For lower values of X2 the 
ratio is greater than one, and (with large errors) appears to decrease at higher values of 
X2- Where the ratio exceeds one, from Eq. this implies d^{x) > vP{x). Furthermore, 
with this data one can map out the difference between dF{x) and vP{x) over a substantial 
region of x. 

The curves in Fig. |1^ are from phenomenological parton distributions [pO|-|5^. The 
upper curves allow (^{x) ^ vP{x), while the lower curves constrain the up and down 
sea quark distributions to be identical. For X2 < 0.2, the ratio of DY cross sections is 
clearly greater than one. In the following subsection we discuss the implications of the 
E866 results, and we examine several theoretical models which might generate large flavor 
symmetry violation in the nucleon sea. 
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Figure 13: The ratio d^{x)/u^\x) extracted from preliminary E866 results of Fig. 12 



4.2 Implications of Large Flavor Symmetry Violation from Drell- 
Yan Experiments 

The experimental Drell-Yan results appear to show a substantial violation of flavor sym- 
metry in the proton sea. The preliminary results of FNAL experiment E866 |]T5[ are 



most definitive in this regard, as they can extract (P'{x)/vP{x) over a substantial range 
of X. In Fig. |13] we plot the ratio (F{x)/vF{x) vs. x which has been extracted from the 
(preliminary) E866 data |jl5|. 

The FSV contribution seen in the Drell-Yan experiments is surprisingly large, as 
it is much larger than can be accommodated by perturbative QCD. Both NLO and 
NNLO QCD calculations have been carried out, and predict very small FSV effects. 



Ross and Sachrajda showed that the flavor symmetry violating contribution which 
arises from higher order QCD evolution is very small (for a comprehensive review of 
theoretical estimates of FSV, see the review by Kumano [^). Consequently, we need 
a non-perturbative mechanism to generate flavor-violating sea quark distributions which 
will reproduce the experimental result. 



Several authors (see the review by Kumano |46] for references) have investigated 
mechanisms for producing a large excess of d over u in the proton. Since there are two 
valence up quarks and one valence down quark in the proton, the Pauli principle should 
make it easier to form a dd pair than a uu pair in the presence of the valence quarks 
(one should be somewhat careful of statements like this: a recent paper by Steffens and 
Thomas suggests that uu pairs may in fact be favored if all antisymmetrization terms 
are considered). In Feynman and Field's early paper on parton distributions |5^, they 
assumed an excess of d quarks in the proton on these grounds. 

Another mechanism for generating additional d quarks in the proton is due to the 
"Sullivan Effect" [^, whereby the virtual photon can couple to a meson created by a 
quark. Thomas |EB| pointed out in 1983 that "pionic" effects would produce an excess 
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Figure 14: Schematic picture of "mesonic" contributions to DIS for a nucleon. a) virtual 
meson contribution; b) recoil baryon contribution. 



of d over u in the proton. The basic mechanism is shown in Fig. [T^. A proton fragments 
into a neutron and a 7r+. If the virtual photon scatters from the down antiquark in the 
Tr"*", it will produce an excess of d over u. 

Several authors have shown that effects given by the "Sullivan" process can produce an 
excess of d quarks in the proton sea. Detailed reviews of this process, and of the literature 
on this subject, are given by Kumano [^], and by Speth and Thomas The original 
papers ||5^ considered the contributions to this process from excess pions in the nucleus. 
Kumano and Londergan |jlT|] calculated a model which included contributions from pions, 
nucleons and A isobars (the A contribution tends to cancel the nucleon-only contribution 
at larger x). This relatively simple model has been expanded by the Adelaide |^ and 
Jiilich 1^ groups to include all of the meson and baryon states normally associated with 
"meson-exchange" models. 

The "mesonic" models look promising, in that both the magnitude of this effect, 
and the x-dependence of the predicted d^/vP difference, are in qualitative agreement 
with experiment. In Fig. |13|, the solid curve is the CTEQ4M phenomenological parton 
distribution |5^. This is quite similar to the 'mesonic' model result of the Jiilich group 
|l58| for the ratio of down to up antiquark distributions. The model of Kumano and 
Londergan which has only vr and A in addition to the nucleon, is similar to the 
CTEQ4M prediction at small x but gives a smaller ratio at large x. Both mesonic 
calculations correctly predict the excess d^{x) > vP{x), and both get the general shape of 
the down antiquark excess as a function of x. At larger x the error bars are sufficiently 
large that detailed comparisons are difficult. Furthermore, the mesonic models are very 
sensitive to small changes in the nNA coupling constant, and to the shapes of the NNir 
and A^Avr form factors. 

There are also other theoretical models which predict an excess of d antiquarks in 
the proton. Eichten, Hinchliffe and Quigg |]59| have investigated the contribution from 



a model in which quarks couple chirally to pions 0. Dong and Liu estimate the 



*We note that this tends to overestimate the asymmetry as it overlooks constraints on the quark 
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contributions from mesons in lattice gauge calculations. They try to separate the contri- 
butions from "cloud" antiquarks and "sea" antiquarks in a lattice calculation. It is not 
possible to make a precise separation on the lattice, but their calculations also suggest 
an excess of down antiquarks relative to up antiquarks. All of this work is summarized in 



Kumano's review article |46]. One additional possibility is that instanton condensates in 
the nucleon pT| , |52| might produce an excess of down sea quarks relative to up sea quarks 
in the proton. 

The Drell-Yan experiments appear to show conclusively a large violation of flavor 
symmetry in the proton sea. However, it is important to note that all these results 
depend on the assumption of parton charge symmetry. If one relaxes this assumption, 
one could in principle reproduce the Drell-Yan results even if flavor symmetry is exact. 
From Eq. let us assume exact flavor symmetry in the proton sea, i.e. 

uP{x) = dF{x) = fix) 

qP{x) ^ (58) 

The parton distributions of Eq. ^ are completely flavor symmetric but not charge sym- 
metric. If we go to the region x^? > we find that Eq. ^ now becomes 

The Drell-Yan experiments could thus be reproduced, even if flavor symmetry was exact, 
with a sufficiently large violation of charge symmetry in the parton distributions. It would 
require an astonishingly large CSV contribution in the nucleon sea to reproduce the E866 
results: this would be a factor 25-50 larger than our estimates in Sect. 3. Alternatively, 
the E866 results could be due to a linear combination of FSV and CSV effects in the 
nucleon sea. 

In the next subsections, we will review the experimental constraints on charge sym- 
metry in parton distributions. In Sect. 5, we will suggest a number of new experiments 
which might provide more stringent tests of parton charge symmetry. 

4.3 The "Charge Ratio:" Comparison of Muon with Neutrino 
Induced Structure Functions 

In Sect. 2.5, we derived the relation between the structure function F2 measured in 
neutrino induced charged current reactions, and the F2 structure function for charged 
lepton DIS, both measured on isoscalar targets. From Eq. |3^, at sufficiently high energies 
the structure functions have the form 



F^^%x) = X 



5 , 1, , ^ , ^^ 4:6d(x)+4:6d(x) + 6u(x) + 6u(x) 
Q{x) - -{s{x) + s{x)) ^ ^ ^ ^ ^ ^ ^ ^ 



18 ' 6 ' ' ' ' 18 



states available in the hadron spectator. 
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i'2 [X) = X 



, 5d{x) + 5d{x) + 5u{x) + 5u{x) 



j=u,d,s,c 

^ _ i(£M±£W) (60) 

In Eq. ^ we have neglected the charmed quark contribution to the structure functions, 
and for the moment we have set R = 0. The function F2^^°{x) is the average of the 
neutrino and antineutrino induced charged current structure functions. 

From Eq. ^ we see that there is a simple relation between the two structure functions, 
in the limit of exact charge sjTumetry. The ratio of the two structure functions in Eq. |6C| , 
when corrected for the strange quark contribution and the factor 5/18 (which reflects the 
fact that the virtual photon couples to the squared charge of the quarks while the weak 
interactions couple to the weak isospin), is defined as the "charge ratio" Rc{x, Q^) or, as 
it is sometimes termed, the "5/18*'^ rule." This quantity should be one, independent of 
X and Q^, in the naive parton model. If we expand the ratio Rc to lowest order in the 
presumably small charge symmetry violating terms, we obtain 



AF2^^°(x)-a;(s(a;) + s(x))/6 



RJx) = 

18 

3 (6u{x) + 5u{x) — 5d{x) — 5d{x) 

^ 1 + ^ ^ ^ (61) 

10Q{x) 

As we pointed out in Section 2.6, the strange quark distributions can be obtained inde- 
pendently by measuring opposite sign dimuon events in neutrino DIS from nuclei. Using 
these strange quark distributions in Eq. |6l|, and comparing the F2 structure functions for 
lept on- induced processes with the F2 structure functions from weak processes mediated 
by W^-exchange, we can in principle measure parton charge symmetry violation and de- 
termine its x dependence. By measuring Rc{x) we can place upper limits on parton CSV 
as a function of x. The longitudinal/transverse ratio R can be included in forming the 
structure functions, and will cancel when the ratio Rc is taken. 

Eq. ^ requires averaging the F2 structure functions for neutrino and antineutrino 
cross sections. If we instead take the ratio using only neutrino-induced reactions, it is 
straightforward to obtain 



R'Jx 



f:I''\x) 



cv I ^fY^''''{x)-x{s{x)+s{x)) /Q 

s{x) - s{x) (45m(x) - 5u{x) - A5d{x) + 5d{x] 



Q{x) 5Q{x) 



(62) 



Eq. ^ differs from Eq. |^ since it has a term proportional to the difference between 
strange and antistrange quark distributions, and also in the relative weighting of the 
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CSV terms which enter. The s — s term is absent if one is able to average neutrino and 
antineutrino cross sections. 

The charge ratio test allows us to place the strongest limits to date on parton CSV. 
There should be no additional QCD corrections to this relation so it should be inde- 
pendent of Q^, provided thai the structure functions are calculated in the so-called "DIS 
scheme." In this scheme, the F2 structure functions are defined so that they have the form 
^2(0;) = xJ^i^'iiQiix) + qi{x)] to all orders, where Cj is the quark charge appropriate for 
either the electromagnetic or weak interactions. For example, the CTEQ4D parton dis- 



tributions ||50|] were determined in the DIS scheme. Despite the robustness of the charge 
ratio test, it also depends on a large number of assumptions and corrections, which must 
be taken into account to obtain limits on CSV terms. Among these corrections are: 

• Relative normalization between leptonic and neutrino cross sections. 

• Corrections due to strange quarks. As outlined in Sect. 2.6, s{x) {s{x)) can be 
independently extracted from the cross section for opposite sign dimuons from 
reactions induced by neutrinos (antineutrinos) 

• Corrections from excess neutrons. Eq. ^ was derived for isoscalar targets. In order 
to obtain reasonable cross sections, neutrino reactions are now measured on iron 
targets. This requires a correction for the excess neutrons in the target. 

• Heavy target corrections. If the leptonic structure functions are obtained from light 
targets and neutrino reactions performed on heavy targets, it is necessary to correct 
the neutrino F2 structure functions for heavy target effects. At low and intermediate 
X, heavy target structure functions are decreased because of shadowing and EMC 
effects, respectively; at very large x Fermi motion effects increases the structure 
functions for heavy targets. 

• Higher twist effects on parton distributions. 

• Heavy quark threshold effects. At sufficiently low energies, heavy quark threshold 
effects will modify structure functions, as we reviewed in Sect. 2.2. 

In Fig. |15| we plot the charge ratio W^, i.e. the ratio of muon F2 structure functions 
measured by Meyers et al. on iron to the value of F2 extracted from the CCFRR 
neutrino measurements [Q. The muon measurements were taken at FNAL with 93 and 



215 GeV muons, using the multimuon spectrometer at FNAL. The CCFRR neutrino 
measurements were made with the FNAL narrow-band neutrino beam. 

In comparing the muon and neutrino measurements, the following corrections were 
made by Meyers et al.. First, the F2 structure functions were modified by including 
the strange quark contribution, determined as described in Section 2.6 [^. Second, 
corrections were made for the excess neutrons in iron. Third, there was a discrepancy 
in the extraction of the F2 structure functions. The muon data were analyzed assuming 
longitudinal/transverse ratio i? = 0, while the neutrino data assumed R = 0.1. Meyers 
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Figure 15: Charge ratio R^ix) of Eq. |6 
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vs. X. Muon measurements of Meyers et al, 
using muons on Fe at FNAL], with CCFRR neutrino measurements on iron, 
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et al. corrected the muon data to make them consistent. The muon data have been 
renormahzed by the factor 1.025. 

Within the errors (two standard deviations), i?^ is consistent with unity, except pos- 
sibly at the largest value of x where R^ix = 0.65) = 0.82 ±0.09. From these experiments, 
the upper limits on the CSV contribution to Rc are generally no better than about 10%, 
and at large values of x the errors are significantly larger. The experimental data is 
consistent with zero charge symmetry violation and certainly rules out any extremely 
large violation of parton charge symmetry. From Eq. ^ and the theoretical calculations 
of parton CSV given in Sect. 3, we expect that the CSV contribution to the charge ratio 
will not exceed 1-2% at any value of x. Consequently, any deviation of the charge ratio 
from unity, at any value of x, would be surprising and very interesting. 

More recent data for both muons and neutrinos allows us to make substantially more 
precise tests of parton charge symmetry. The NMC group Q measured the F2 structure 
function for muon interactions on deuterium at energy i?^ = 90 and 280 GeV. The CCFR 
group has extracted the F2 structure function for neutrino and antineutrino interac- 
tions on iron using the Quadrupole Triplet Beam at FNAL. The CCFR measurements 
provide the most copious sample of neutrino events, and allow the most precise limits 
on parton CSV. In Fig. |16| we plot the charge ratio Rc of Eq. ^ vs. x. The circles are 
the NMC/CCFR ratio. The open triangles are the BCDMS/CCFR charge ratio, where 
BCDMS represents the muon scattering results of the BCDMS group on deuterium 
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Figure 16: Charge ratio R^ix) of Eq. |6^ vs. CCFR neutrino measurements on iron 
with the FNAL wide-band neutrino beam, Ref. Circles: ^ + D measurements of 

the NMC group, Ref. open triangles, muon measurements from the BCDMS group 
on deuterium, Ref. and carbon, Ref. solid triangles: SLAC electron scattering 



data, Ref. [24,68 



and carbon |6^. The solid triangles are the SLAC/CCFR charge ratio, where SLAC 



denotes electron scattering results of the SLAC group 



The charge ratio has been calculated by C. Boros |69|. The results differ somewhat 



from those produced by Seligman et al. in their calculation of the charge ratio [|65| , |f0[ . In 



comparing the data sets, Boros takes only those points with the same x value and sums 
over overlapping values, while Seligman interpolates between measured values of the 
structure functions. In addition, in Fig. ^ there is no correction for strange quarks. 

In the region 0.1 < x < 0.4, the charge ratio test is consistent with unity, and the 
data gives an upper limit to CSV effects in the charge ratio at about the 3% level. 
For larger values of x the upper limit on CSV effects is more consistent with the 5- 
10% level, due mainly to the poorer statistics and, as we will see, on the large Fermi 
motion corrections needed for the heavy target at large x. Both the new muon and 
neutrino data are more precise than the older measurements. In addition, the more recent 
phenomenological parton distributions are better determined. Relative normalizations of 
lepton and neutrino cross sections appear to be well understood. All data is analyzed 
with consistent assumptions about the longitudinal/transverse ratio R. Heavy quark 
threshold effects should also be under control. 

Probably the most significant correction is the heavy target correction, necessary 
because we are comparing muon data on deuterium, where the correction is presumably 
very small, to neutrino data on iron. In Fig. we show the same charge ratio as in 
Fig. |1^ for the NMC-CCFR comparison, but here we explicitly show the heavy target 
corrections. The open triangles show the ratio without heavy target corrections, and the 
solid circles show the ratio after applying these corrections. The solid line is the iron 
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Figure 17: Effects of heavy target correction on charge ratio of Fig. |T6|. Triangles: NMC- 
CCFR data without heavy target correction; circles: NMC-CCFR data including heavy 
target correction (solid curve). 



target correction factor as a function of x. 

After including the heavy target correction, there appears to be a significant deviation 
of the charge ratio from one at the smallest values x < 0.1; the discrepancy approaches 
15% at the smallest values of x, with the electromagnetic structure functions being smaller 
than the neutrino ones. Several suggestions have been made to explain this discrepancy. 
We summarize the explanations listed by Seligman [Q. First, the lack of agreement 
could result from difficulties in analyzing the low-x neutrino events. This will be tested 
with the next generation FNAL neutrino experiment E815. Second, it is conceivable 
that the disagreement arises from effects at small which differ between leptonic and 



neutrino induced reactions fTIl]. However, these effects appear to be quite small for 
> 1 GeV2 0. 

The discrepancy increases monotonically at small x, where the strange quark effects 
are largest. One intriguing possibility is that strange quark effects might account for 
all of the apparent discrepancy. In this case it is possible that the present phenomeno- 
logical analysis of both strange quark and antiquark distributions need to be modified 
substantially, as has recently been argued by Brodsky and Ma [^. In any case, the 
recent NMC-CCFR comparison allows us to put rather tight upper limits on parton CSV 
contributions, and focuses our attention on the low-x region where there is currently a 
discrepancy between the F2 structure functions extracted from the two reactions. 
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4.4 Comparison of Neutrino and Antineutrino Cross Sections 
on Isoscalar Targets 

On an isoscalar target, the differential cross sections for charged current reactions induced 
by neutrinos or antineutrinos can be written in the general form 



dx 
dx 



G^MnE 
n 

TT 
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+ 



1 



2 ' Q^i+R)]^^"'''^^^Q')+r^^"''"(^^Q') 
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2 6(1 + 

In Eq. ^ the quantity R is the longitudinal/transverse ratio, 

^2, CTL F2(x,g2) _2xFi(x,g2 



(63) 
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(64) 



From Eq. ^ we see that the structure functions F2 and F3 per nucleon for an isoscalar 
target can be written as 
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u^ix) + dP{x) + vP{x) + dP{x) + 2s{x) 
vP{x) + dP{x) + vF{x) + dF{x) + 2s{x) 
uP{x) + dP{x) - vP{x) - dP{x) + 2s{x) 
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— 6u{x) — 6d{x) 

— 6d{x) — 6u{x) 

— 6u{x) + 6d{x] 

— 6d{x) + 6u{x] 



(65) 



In Eq. ^ we assume that the momentum transfers are sufficiently high that threshold 
effects can be neglected. In this equation, we have neglected the contribution from 
charmed quarks in the nucleon, and for the moment we have set i? = 0. In this limit, 
the F2 structure functions from neutrinos and antineutrinos are identical except for CSV 
contributions. In addition, the F2 and F3 structure functions are identical except that 
the antiquark contributions have different signs. Consequently, if we go to large x where 
the sea quark contributions become small with respect to the valence quark terms, then 
both F2 and F3 structure functions for both neutrinos and antineutrinos should become 
equal. From Eq. |6^ we see that F2 and F3 will add together in the neutrino cross section, 
but will cancel in the antineutrino cross section. 

We thus define the ratio of antineutrino to neutrino charged current cross sections on 
an isoscalar target, 

^ da^^^(x)/dx , , 

^ ' da''^"{x)/dx ^ ' 

We will focus on this relation at reasonably large values of x > 0.3. For these values of 
X the sea quark distribution will be small relative to the valence quark distributions. In 
this region we can expand the ratio of Eq. to lowest order in small quantities, and we 
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Figure 18: The quantity 1 - Sr'^/^ of Eq. |66|, vs. X for neutrino and antineutrino data on 
BEBC D bubble chamber. Data from WA25 experiment, Ref. W^. 
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In Eq. we treat as a constant (we use the value averaged over x). If i? = 0, the 
ratio r'^/^(x) is predicted to approach the value 1/3 at large values of x; for R = 0.1, the 
ratio should approach 21/65. x-dependent deviations from this constant value will arise 
from either sea quark or CSV contributions. 

At present, this ratio can be obtained from experimental data on deuterium and iron 
targets. In Fig. 18 we plot 1 — 3r'^/'^ vs. x for the WA25 data ||7^. This consists of neutrino 
and antineutrino cross sections taken in the BEBC D bubble chamber at CERN. The 
experimental points plotted have an average momentum transfer = 11 GeV^. They 
were analyzed assuming R = 0. The dotted curve in Fig. |l^ is the total contribution from 
both sea quarks and our model predictions from valence quark CSV (see Section 3), while 
the dashed curve is the model contribution from CSV alone. The model prediction is 
small in absolute value, and small relative to the sea quark contribution, except at large 
values X > 0.6, where the CSV contribution is predicted to dominate. In the quantity 
1 — 3r, the sea quarks are weighted with a factor 8/3 relative to the CSV contribution. In 
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Figure 19: The quantity 1 - 65r^/^/21, vs. X for CCFRR neutrino and antineutrino data 



on iron, Ref. IH^ . 



this ratio, we predict that the sea quarks will give a negative contribution to the quantity 
1 — 3r, while the CSV contribution is predicted to be positive. 

The experimental error bars range from about 20% at small x to 50% at large x. The 
errors are substantially larger than the theoretical CSV contribution, for all measured 
values of x. If it were possible to obtain precise neutrino and antineutrino data on isoscalar 
targets such as D, for large values of x, the quantity could in principle give a 

sensitive upper limit on parton charge symmetry violation. Unfortunately, measurement 
of these cross sections is notoriously difficult. We know of no current plans for precision 
measurements of charged current cross sections for neutrino and antineutrino beams on 
isoscalar targets. From Fig. 0, we see that the experimental error bars would have to be 
at least an order of magnitude smaller than their current values to reach our predicted 
CSV signal. Further, it is doubtful that the neutrino results could attain the limits on 
CSV already reached by the "charge ratio" comparison between F2 measurements from 
muons and neutrinos; this was reviewed in Sect. 4.3. 

This is unfortunate, since the neutrino comparisons have fewer implicit assumptions 
than the charge ratio - if accurate structure functions were available for neutrino and 
antineutrino bombardment of isoscalar targets, one would be comparing data taken in 
the same experiments, one would not have to make corrections for excess neutrons, nor 
would one have to correct for the strange quark distributions, extracted as described in 
Sect. 2.6. 

In Fig. PI we plot 1 - Q^r"/^ /2l vs. X for the CCFRR data, neutrino and antineutrino 
cross sections on iron, taken at FNAL |Q. The plotted points correpond to an average 
momentum transfer = 12.69 GeV^. They were analyzed assuming R = 0.1. The 
dotted curve in Fig. |T9| is the total contribution from both sea quarks and our model 
predictions from valence quark CSV, while the dashed curve is the model contribution 
from CSV alone. The CCFRR data is similar to the WA25 results shown in Fig. 113, in 
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that the experimental error bars are much larger than the theoretical CSV contribution, 
and the data is consistent with no charge symmetry violation. 

The experimental error bars range from about 20% at small x to 100% at large x. 
The errors are again substantially larger than the theoretical CSV contribution, for all 
measured values of x. Since iron is not an isoscalar target, it is necessary to make 
corrections from the excess neutrons in iron. These corrections have been taken into 
account in Figs. and |19|. There is more recent experimental data for neutrino and 
antineutrino cross sections from the CCFR collaboration We are unaware of a 

systematic study of the quantity r'^^'^{x) by this group. 



Proposed New Experimental Tests of Parton Charge 
Symmetry 



In the preceding Section, we reviewed existing experiments and showed the limits 
they placed on charge symmetry and flavor symmetry violation in parton distributions. 
The latest Drell-Yan data for protons on proton and deuteron targets appear to show 
clear evidence of substantial flavor symmetry violation in the proton sea. However, it 
is conceivable (although unlikely) that this result could also be due to charge symmetry 
violation in parton distributions, or by some combination of flavor symmetry and charge 
symmetry violation. In this Section we propose a series of experiments, all of which are 
chosen specifically to set limits on CSV contributions to parton distributions. 

5.1 Test of Weak Current Relation F^^^"{x) = 

In Sect. 2.5, we reviewed the high-energy limiting form to the electroweak structure 
functions. We showed that, at sufficiently high energies, the charge- changing structure 
functions on an isoscalar target are equal except for contributions from valence quark 
CSV, and possible strange or charmed quark terms, i.e. 



as shown in Eq. ^ 

At the enormous values of that can be probed at HERA, weak interaction processes 
such as e^p u^X are not completely negligible with respect to the electromagnetic 
process e~p — > e~X. Furthermore, we are well above heavy quark production thresholds, 
so threshold effects on the CKM matrix elements and issues like "slow rescaling" which 
are critical at lower energies, are unimportant in this regime. The (e~, Ue) reaction picks 
out the positively charged partons in the target, so that if one could accelerate deuterons 
in the HERA ring, the structure functions for this reaction would have the form 

F^'^ix) = lu^ix) + dPix) + + r{x) + 2s{x) + 2c(x)l . (69) 
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As in Sect. 2, we denote the (e~,i/e) reaction by the charge of the virtual W absorbed 
by the target. We neglect differences in strange and charm quark distributions between 
proton and neutron. On the other hand, if we have a positron beam the (e"'",z/e) deep 
inelastic reaction measures only the negatively charged partons, so for a deuteron target 
this would have the form: 

F^^^ix) = [d^ix) + uP{x) + c/"(x) + + 2s{x) + 2c{x)] . (70) 

Taking the difference of the e"*" and e~ charge- changing weak interaction cross sections 
one therefore obtains Eq. The difference between the structure functions 
and Fj^ ^ has been studied recently by Londergan, Braendler and Thomas [Q. We 
summarize their results here. 

To indicate the size of expected differences in the e^D charge-changing cross sections, 
we construct the ratio 



Rwix) 



2 [Fi^^^ix) - Fl 



F^^D{x) + F^- 



X] 



Sdv{x) — 6uv{x) + 2 [s{x) — s{x)] 



J2j=p,n u^{x) + ui{x) + di{x) + di{x) + 2{s{x) + s{x)) 
= Rcsvix) + Rsix) (71) 

Direct comparison of the Fi structure functions for charge-changing weak interactions 
on an isoscalar target is a strong test of charge symmetry in parton distributions. As 
we have shown, the only source of difference is either a difference between strange quark 
and antistrange quark distributions, or charge symmetry violating components of the 
valence quark distributions. Note that Eq. |7l] is unchanged if we use the F2 structure 
functions rather than the Fi structure functions: the additional factors proportional to 
the longitudinal/transverse ratio R cancel in forming the ratio. Eq. is also true for any 
isoscalar nuclear target, if we replace the nucleon parton distributions by their nuclear 
counterparts. 



In order to illustrate the size and shape of the effect expected we plot in Fig. 20a the 



theoretical CSV contribution, Rcsv{x) from Eq. The dashed curve is calculated for 
Q"^ = 100 GeV^, the dot-dashed curve for Q"^ = 400 GeV^, and the dash-triple dot curve 
for Q2 = 10,000 GeV^ The CSV quantity Rcsvix) is predicted to be greater than 0.02 
provided x > 0.4. The valence quark CSV terms are taken from the work of Rodionov et 
ai, Ref. [|^, as discussed in Sect. 3.1. For values x > 0.1 we predict that 6dv{x) will be 



positive and 6uv{x) negative, so their effects should add. As we mentioned previously, 
estimates of valence quark CSV have little model dependence, so we would expect to see 
differences of this magnitude and sign at large Bjorken x. We predict several percent 
effects at the largest values of x, if the structure functions could be determined in this 
region. 

The term Rg of Eq. |7l| is shown in Fig. pO|b. The term Rs{x) is proportional to the 
difference between strange quark and antiquark distributions. There has been quite a lot 



46 



iOO GeV ' 
400 GeV ' 
10000 GeV ' 



0.0 0.2 0.4 O.f 



■K \\ 
I 


. 1 00 GeV ' 


\'; 
V-; 


400 GeV ' 

10000 GeV ' 


\ 
\ 




\ 
\ 




\ 




\ 
'■\ 
w 

\ 




\ 


\ 
\ 

. . 1 . < . 1 . . 





0.0 0.2 0.4 0.6 



Figure 20: a) Theoretical estimates of CSV contribution Rcsv to the ratio Rw{.x) of Eq. 



|7T| vs. X, for various values of . 
Eq. |7l|. From Ref. WA\. 



b) Theoretical estimates of s — s difference, Rs{x) of 



of interest recently p4^36| in the possibility that s{x) — s{x) might be non-zero. The 
"mesonic models," which have had success in reproducing the experimental values for 
d{x)/u{x), as discussed in Sect. 4.1 and 4.2, naturally give rise to differences in these 
distributions. In such models, the s arises from virtual kaon production, while the s 
comes from the residual strange baryon (A or S); this was first pointed out by Signal 



and Thomas |3^. There was some suggestion of experimental support for the idea |28 
based on LO analysis of the CCFR neutrino data, but subsequent NLO analysis of the 
same data ||30| saw no difference between strange and antistrange quark distributions, 
within the experimental errors. We reviewed the current experimental situation regarding 
strange/antistrange quark distributions in Sect. 2.6. 

For an estimate of the difference between strange and antistrange quarks, we have 
used the mesonic model calculation of Melnitchouk and Malheiro For this estimate 



we took the (poorly determined) NKH vertex function to be a monopole form factor 



of mass 1 GeV, the largest value consistent with the latest CCFR data |2^, Fig. 
?? shows the ratio Rs of Eq. ^ from this model. The order of magnitude of the s — s 
difference is comparable to that arising from CSV. As the two effects have the opposite 
sign, we predict some cancellation between the two contributions. However, the predicted 
shapes are completely different, and one should be able to separate the two contributions 
on the basis of the measured x-dependence of Rw- 

We emphasize that even the sign of the s — s difference is not well determined, so the 
theoretical "error bars" associated with the curves of Fig. ?? are large. As we mentioned 
in Sect. 2.6, the quantity s — s can be independently extracted by measuring opposite 



47 



sign dimuons arising from scattering of neutrinos or antineutrinos from nuclei. 

The structure functions F^^^ and ^ should be nearly identical at all x. If 
they are not, this would be quite surprising. Any deviation from zero, at any value 
of X, would he extremely interesting, whether its origin lies in parton CSV or intrinsic 
strangeness. This is a strong test of the validity of charge symmetry in parton distri- 
butions. Since the ratios require comparison of charge- changing reactions induced by 
electrons and positrons, it is important to have very accurate calibration of the relative 
reaction rates. Detector efficiencies should not be a major problem as the signal involves 
prominent jets on the hadron side and very large missing energy and momentum on the 
lepton side. Clearly it will be important to determine experimentally whether or not 
pW+D _ pw D non-zero. The interpretation in terms of CSV, s 7^ s or possibly both, 
can then be pursued in detail. 

5.2 Drell-Yan Processes Initiated by Charged Pions on Isoscalar 
Targets 

In Sect. 4.1 and 4.2, we showed the dramatic results obtained by comparing pp and 
pD Drell-Yan [DY] processes. The preliminary results from the E866 experiment at 



FNAL [|T5| appear to show substantial flavor symmetry violation in the proton sea. We 
can also use DY processes as a specific test for charge symmetry violation in parton 
distributions. For this we want to differentiate between, say, up quarks in the proton and 
down quarks in the neutron. We will show that this could be accomplished by comparing 
DY processes induced by charged pions on isoscalar targets. 

The crucial element here is that at large momentum fraction x, the nucleon distri- 
bution is dominated by its three valence quarks, while at similar large x the pion is 
predominantly a valence quark-antiquark pair. If one uses beams of pions, and concen- 
trates on the region where Bjorken x of the target quarks is reasonably large, then the 
annihilating quarks will predominantly come from the nucleon and the antiquarks from 
the pion. Comparison of Drell-Yan processes induced by 7r+ and 7r~ in this kinematic 
region provides in principle a sensitive method for comparing d and u quark distributions 
in the nucleon, since the 7r+ contains a valence d (and will annihilate a d quark in the 
nucleon) and vr" a valence u (and will annihilate a nucleon u quark). 

As an example, we consider reactions on the deuteron, although our results will be 
true for any isoscalar nuclear target. Consider the DY process for a 7r+ on a deuteron. In 



Fig. |TT]b, we showed a schematic picture of the dominant process (in a kinematic region 
dominated by valence quarks for the meson and nucleon), for charged pion DY processes 
on a proton. Provided that x, > 0.3, to minimize the contribution from sea quarks, 
the dominant process will be the annihilation of a c? in the vr"*", with momentum fraction 
Xtt, with a down quark in the deuteron with momentum fraction x. Neglecting for the 
moment sea quark effects, the tx^ — D (vr^D) DY cross sections will be proportional to: 

(T^+^(x,x^) ~ ]-{d'P{x) + d''{x))'d!' (x„) , 
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a^''j,{x,x^) ~ UuP{x)+u^{x))lf'{x^) . (72) 



If we construct the ratio, R^^{x,Xt,): 

-"^ttD l-^' -^ttJ ~ /A DY ( \ I DY ( W Ir^i V^) 



{4a^rDi^,x^) + a^^^{x,x^)) /2 ' 



we expect the ratio R^^ to be quite sensitive to charge symmetry violating (CSV) terms 
in the nucleon valence parton distributions. 

From our model calculations of Section 3, the CSV contributions should be at most 
a few percent. Consequently we must include sea quark contributions for both nucleon 
and pion in defining this ratio. There will also be a contribution from charge symmetry 
violating effects in the pion parton distributions. Including these contributions, the Drell- 
Yan ratio for pions can be written: 



+ 



/ 6d{x) — 6u{x) \ 
[uZ{x) + d^.{x) ) 

15 [27r,(x^)<(x) + TTsjx^) {uPjx) + dPjx))] _ 

4(i^(x7r) [uZ{x) + d?,{x)] 
<^(x) + i?f^(x,x.), (74) 



Eq. ^ is valid at sufficiently large x and x^, where sea quark probabilities are small 
relative to valence quarks. We have expanded it to lowest order in both sea quark and 
CSV terms 0. 

The nucleon CSV term i?^(x) in Eq. (^) is a function only of x. It is not necessary 
to know absolute fluxes of charged pions to obtain an accurate value for i?^^. The yield 
of J/ip^s from TT"*" — Z) and t[~ — D should be identical to within 1%, so this can be used 
to normalize the relative fluxes. Because R^^{x, x^) is a ratio of cross sections, a number 
of systematic errors should cancel. In particular, Eq. (|74D is not sensitive to differences 
between the parton distributions in the free nucleon and those in the deuteron |^,|7S|-|S^. 



Provided that both the neutron and proton parton distributions are modified in the same 
way, then the ratio in Eq. ([7^ ) will be unchanged. 

In Fig. ^ we show the nucleon CSV contribution, R^^{x), using the bag model 
predictions for CSV, evolved to = 10 GeV^. As the main uncertainty in our calculation 
is the mean diquark mass, the results are shown for several values of this parameter. In 
the region 0.4 < x < 0.7, we predict R^ will always be positive, with a maximum value 
of about 1.7%. 

We predict that the contributions from pion and nucleon CSV will all be the same 
sign and will add constructively. 



5- — + 



in Eq. |74| we neglected a pion CSV term Sd (xtt) = d (x^) — (xtt). This term was estimated 



using a Nambu-Jona Lasinio (NJL) ||75| model employed recently by Toki and collaborators 76, 7' 
which predicted a very small pion CSV effect 0]. 
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Figure 21: Theoretical estimate of nucleon CSV term R^j^ of Eq. From Ref. 



Since the Drell-Yan ratios arising from CSV are expected to be very small, even small 
contributions from sea quarks could make a substantial effect. The dominant contribution 
will arise from interference between one sea quark and one valence quark. We concluded 
that for sufficiently large x, e.g. x > 0.4, one should be able to separate the CSV "signal" 
from the sea- valence interference |^. Unlike the CSV contributions of Eq. (|7^), the sea- 
valence interference term R^^{x,Xt,) does not separate, so one could exploit the very 
different dependence on x and x-,, of the background and CSV terms. We conclude that 
the CSV terms could be extracted even in the presence of a sea- valence "background" . 

Despite our prediction in Sect. 3 that the fractional "minority quark" CSV term, 
6d{x)/d^{x), should be between 3 and 7% (c.f. Fig. |10D, the nucleonic CSV ratio R^d, 
shown in Fig. |2T|, is predicted to be more like 1-2%. This is because 6d in Eq. ( [TID is 
divided by + dP, and at large x d^{x) << u^{x). A much larger ratio could in principle 
be obtained by comparing the vr"^ — p and "vr" — n" Drell-Yan processes through the ratio 



4cTf/ (x, x^) + erf J (x, x^) 



a. 



DY 
-D 



4^f+;(^, X.) - <%{x, x^) + af4(x, X.) /2 



(75) 



In principle, the advantage of this measurement is that it isolates the minority quark 
CSV term - in fact, the dominant term in Eq. |75| is the term 6d^{x)/d^{x) so we expect 
CSV effects at the 3-7 % level. 

We conclude, however, that this quantity is unlikely to provide unambiguous infor- 
mation regarding parton charge symmetry violation. First, to form the ratio in Eq. |75| 
one must know the relative normalization of DY cross sections on protons and deuterons. 
This should be feasible by bombarding both hydrogen and deuterium targets simulta- 
neously with charged pion beams. In order to extract the minority quark CSV term, 
it is necessary to know the precise relation between deuteron parton distributions and 
those for free nucleons. If we include "EMC" changes in the deuteron structure func- 
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tions relative to free proton and neutron distributions, we find that 2-3% changes in the 
parton distributions can produce 10-30% changes in the ratios of Eq. In view of the 
sensitivity of this ratio to the EMC term, we conclude that information regarding CSV 
effects could not be extracted from comparing DY cross sections for tt~p with those for 
"vr"'' — n" unless the Fermi motion and binding corrections for the deuteron were known 
to great accuracy. 

Comparing the Drell-Yan yield for vr^ and tt~ on deuterons may provide a means to 
extract the charge symmetry violating [CSV] part of the nucleon parton distribution. As 
the X and values of interest for the proposed measurements are large (x > 0.5), a 
beam of 40-50 GeV pions will produce sufficiently massive dilepton pairs that the Drell- 
Yan mechanism is applicable. A flux of more than 10^ pions/sec. is desirable, so these 
experiments might be feasible when the new FNAL Main Injector becomes operable. 



5.3 Charged Pion Leptoproduction from Isoscalar Targets 



In a recent paper [^], we pointed out that semi- inclusive pion production, from lepton 
DIS on nuclear targets, could also be a sensitive probe of CSV effects in the valence 
parton distributions for the nucleon. 

In the quark/parton model, the semi-inclusive production of hadrons in deep inelastic 
lepton scattering from a nucleon is given by 

~ WTZK- ("6) 



The quantity A^^'^ in Eq. (|76D, the yield of hadron h per scattering from nucleon A^, 
has the form A^^'* = ^^e^qf {x)D^{z), where D^{z) is the fragmentation function for 
a quark of flavor i into hadron h, which depends on the quark longitudinal momentum 
fraction z = Eh/iy, where Eh and u are the energy of the hadron and the virtual photon 
respectively. 

For pion electroproduction on an isoscalar target, (such as the deuteron) charge sym- 
metry relates the "favored" production of charged pions from valence quarks, by 

Nf^l'ix^z) = ANf:^-ix,z). (77) 

In Eq. 77, Nj?^'^ {x, z) represents the yield of 7r~^ per scattering from the deuteron, via 
the "favored" mode of production, e.g., for tc'^ (tt") production, the "favored" mode of 
charged pion production is from the target up (down) quarks. Since the semi-inclusive 
reactions are proportional to the square of the quark charge, there is a relative weighting 
of 4 for 7r+ production. 



Deviations from Eq. will arise from effects due to sea quarks, CSV effects in the 
parton distributions, and contributions from the "unfavored" fragmentation functions. 
The HERMES collaboration at HERA is currently taking experimental data on 



semi-inclusive pion production from hydrogen and deuterium f\ 

^Our description of fragmentation is correct only in the high energy hmit, where hadron production 
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Assuming charge conjugation invariance and charge symmetry for the fragmentation 
functions allows us to write the yield for leptoproduction of a charged pion from a proton 

as 



ATP" (x,z) 



The fragmentation functions have been extracted by the EMC group [p^, p4 
being independently measured by the HERMES collaboration 
We proposed measuring the quantity R^{x,z), defined by 



(7J 



and are 



R^ix^z) 



l-A(^) 
1 + A(^) 

5A{z) 



AN 



Dtt- 



{x,z) -N^''^{x,z) 



ND^+{x,z) -N^^- 
4 [5d{x) - 5u{x)] 



X, z) 



+ 



1 + A(2) 3 [ul{x) + dl{x)\ 

5 (vP{x) + dP{x)) + A,(z) [s{x) + s{x)] /(I + A(2)) 
[mv(x) + (iv(x)] 



(79) 



In Eq. (|79D, we expand to first order in "small" quantities. These are: the CSV 
nucleon terms, 5d{x) and and the sea quark distributions (Eq. (^) is only valid at 

large x where the ratio of sea/valence quark distributions is small). We have neglected 
the CSV part of the fragmentation function; in Ref. ||81| we estimated that this term 
would be quite small. 

The quantity R^{x, z) in Eq. ^ separates into three pieces. The first piece depends 
only on z. This term is a function only of the experimentally measured quantity A(z), the 
unfavored/favored ratio of fragmentation functions. It decreases roughly monotonically 
as z increases. The second term depends only on x, and is proportional to the nucleon 
CSV fraction (relative to the valence quark distributions). The final term in Eq. |79| 
depends on both x and z. It is proportional to the sea quark contributions, so it becomes 
progressively less important at large x. It also contains a term which is proportional to 
the strange/favored ratio of fragmentation functions /S.s{z): we estimate that this term 
is always negligible. 

Experimentally, one needs to measure accurately the x-dependence of R^{x,z) for 
fixed z; in this case the 2;-dependent term will be large (of order one) and constant. The 
sea quark contribution will be large at small x, but should fall off monotonically and 
rapidly with x. So, at sufficiently large x, the sea quark contribution will be negligible 
relative to the CSV term. One then has to extract the small, x-dependent term in Eq. 
from the large term independent of x. As a general rule, the larger the values of 



is so copious that the leading quark fragmentation and target fragmentation completely decouple. It is 
likely that at HERMES energies, sufficiently few hadrons are produced that this picture is inaccurate. 
Monte Carlo simulations of data at these energies could reveal the breakdown of this naive factorization 
picture, in which case the arguments presented here would be applicable only at higher energies 
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Figure 22: x-dependent contributions to charged pion leptoproduction, from Ref. [Bl 



X and z at which data can be taken, the larger the CSV term will be relative to the 
z-dependent term. 

In Fig. 22 we plot our predictions for the x-dependent terms in R^{x, z), at = 10 
GeV^. The long dashed curve is the contribution from nonstrange sea quarks to R^^^{x, z). 



This depends only on x, and is calculated using the CTEQ3M parton distributions 
The short dashed curve is our prediction for the parton charge symmetry violating term, 
R^svi^)i this uses the CTEQ3M parton distributions, plus the bag model prediction for 
valence quark CSV from Londergan et al. 0], as discussed in Sect. 3. The dot-dashed 
curve is our estimate of the strange quark contribution. The solid curve is the sum of 
the three terms. 

For X ~ 0.5, the CSV term is as large as the sea quark contribution, and with 
increasing x (e.g., for x > 0.55), the CSV term dominates the x dependent terms. We 
predict the maximum CSV contribution will be of order 0.02 — 0.04. The x dependent 
contribution shown in Fig. ^ will sit on a large and constant z dependent term. This 
term is predicted to vary between 1.5 and 1, as z goes from 0.4 to 0.8 |ST]. So the CSV 
term is expected to be between 1-4% of the ^-dependent term. 

In the HERMES experiment at HERA, the goal is to make precision measurements 
of the spin structure functions, so the prospect for obtaining very accurate spin-averaged 
charged pion leptoproduction data is excellent. Only data from deuterium targets is 
required; efficient detection of both signs of charged pions is important, but absolute 
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yields are not required as overall normalizations cancel out in the ratio of Eq. ([75|). 



6 Charge Symmetry and Flavor Symmetry Contri- 
butions to Sum Rules 



Sum rules can provide extremely useful information on parton distributions. For 
example, from the quark model the integrals of up and down valence quark distributions 
obey the quark number sum rules, given by the normalization conditions on the quark 
distributions, see Eq. ^ in Sect. 2. 



can 



Two sum rules, the Adler sum rule |^6|] and Gross-Llewellyn Smith sum rule ||87|| , 
be directly related to linear combinations of quark normalization integrals. The Adler 
sum rule is obtained by integrating the difference between the F2 structure functions 
for charged- current interactions of antineutrinos on protons, and that for neutrinos on 
protons. The Gross-Llewellyn Smith sum rule is obtained by summing the F3 structure 
functions for charged current interactions of neutrinos, and antineutrinos, on a proton 
target (the same result is obtained for a neutron target). 

The Gottfried sum rule (GSR) 0] is obtained by integrating the difference between the 
F2 structure functions for (photon mediated) neutral current interactions on protons and 
neutrons. Unlike the Adler or Gross-Llewellyn Smith sum rules, the "naive" Gottfried 
sum rule expectation Sq = 1/3 is obtained only if we assume both charge symmetry 
for parton distributions, and what is frequently called "SU(2) flavor symmetry" in the 
proton sea. That is, in addition to charge symmetry we assume that 'U^(x) = (^{x). 

There has been much recent interest in the Gottfried sum rule, sparked by the rather 
precise measurements of the NMC group These show rather conclusively that Sq 
is substantially less than 1/3. If charge symmetry is valid, this provides information 
on SU(2) flavor symmetry violation [FSV] in the proton sea. In fact, the level of FSV 
needed to agree with the NMC result is surprisingly large. However, as we will see, the 
experimental GSR measurements are actually sensitive to a combination of FSV and 
CSV in the nucleon sea. 

As we have argued in this paper, we would expect CSV effects in parton distributions 
to be no greater than about one percent for sea quark distributions, and for the "majority" 
valence quark distribution. However, we expect the "minority" valence quark distribution 
to exhibit CSV effects of several percent at large x. As we discussed in Sect. 4, the current 
experimental upper limits on CSV are of the order of a few percent for 0.1 < x < 0.4; 
upper limits for CSV are no better than 10% for larger values of x, and for x < 0.1 the best 
experiments to date suggest a violation of parton charge symmetry. It would therefore be 
useful to construct sum rules which could in principle distinguish between CSV and FSV 
effects. In this Section, we will review the Adler, Gross-Llewellyn Smith and Gottfried 
sum rules, clarifying the various assumptions implicit in their derivation (particularly the 
role of charge symmetry in the sum rules). Next, we will propose sum rules which can 
clearly differentiate between charge symmetry violation and flavor symmetry violation in 
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nucleoli sea quark distributions. Finally, we will present existing deep inelastic structure 
functions in the context of these new sum rules, to examine the degree to which existing 
data can constrain limits on charge symmetry and/or flavor symmetry violating effects. 
The discussion of new sum rules follows rather closely the prior theoretical work of 



Ma [jT6|,|T^]. The sum rules we introduce are linear combinations of those proposed by Ma. 
He also pointed out the potential confusion in the literature on the question of charge 
symmetry. 

We first review existing sum rules (Gottfried, Gross-Llewellyn Smith, Adler) without 
making the usual assumptions of charge symmetry in quark distributions. The NMC 
measurements of the Gottfried Sum Rule are seen as strongly suggesting large flavor 
symmetry violation [FSV] in the proton sea quark distributions. The Drell-Yan measure- 
ments carried out by the NA51 and E866 [|l^ groups are regarded as more or less 
definitive proof of large FSV effects. We point out that all three results could in principle 
be explained by large charge symmetry violation in the nucleon sea quark parton distri- 
butions (alternatively, a linear combination of FSV and CSV effects could be responsible 
for these results). 

In Sect. 7, we introduce "flavor symmetry" and "charge symmetry" sum rules, and 
discuss how they can separate CSV and FSV effects in nucleon sea quark distributions. 
We discuss what (if any) experimental limits on CSV and FSV can be drawn from existing 
deep inelastic neutrino scattering experiments. 



6.1 Review of Gottfried Sum Rule 

Here we review the Gottfried Sum Rule 0. We go through this in considerable detail 
so that the underlying assumptions in its derivation are clear throughout. For a compre- 
hensive review of both experimental and theoretical aspects of the Gottfried Sum Rule, 
and the related question of the flavor symmetry of the proton sea, see the recent work 
by Kumano . 



The Gottfried Sum Rule is given by 



Sg ^ t - [F^%x)-Fr{x)] . (80) 



X 



Because the F2 structure function from electron and muon deep inelastic scattering de- 
pends on the squared charge of the quarks, the up quark distributions are weighted by a 
factor of four relative to the down quark distributions. 

Rewriting the structure functions in terms of quark distributions, we obtain the result 

/' ^ [F^^ix] - F^ix)] = - C dx \Au'\x) +Aunx) - AvJ\x) - Au\x) 
Jo X 9 Jo 

+ (F{x) + dP{x) - (r\x) - J"(x)] . (81) 

In obtaining Eq. (|8TD, we assume the strange quark contributions for neutron and proton 
cancel. We can invoke the "strong" assumption that the strange parton distributions for 
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proton and neutron are identical at each value of x, i.e. s^{x) = (and similarly for 

the antiquark distributions); alternatively, we can assume the "weak" condition that the 
parton distributions need not be identical at all x, but that the integrals over x of the 
appropriate parton distributions are identical for proton and neutron. There is no QCD 
modification of this sum rule. 

Introducing valence quark distributions as in Section 2, i.e. u^{x) = u{x) — u{x), we 
obtain 

r-l dx 1 rl 



Jo X 



r(x) 



Frix)] 



1 

9 Jo 



dx [4<(x) +8mP(x) -4<(x) 



; m"(x) + <(x) + 2 dP{x) - <(x) - 2 c/"(x) 



f82) 



We now invoke the valence quark normalization conditions, defined in Eq. and we 
obtain the Gottfried Sum Rule, 



5', 



G 



dx 



Frix) - Frix) 



3 + 9/0 



Au^ix) + d^ix) - 4m" (x) - d'\x) 



^3) 



If we make the additional (and customary) "strong" assumption of charge symmetry. 



e.g. 



d"(x) 
u'\x) 

rix) 



uPix) 
d^ix) 
uPix) 
d^ix) 



U) 



then we obtain the "normal" formulation of the Gottfried Sum Rule [GSR], 



S, 
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u^ix) - d^ix) 
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Note that the charge symmetric Gottfried Sum Rule, Eq. does not require the "strong" 
assumption of charge symmetry: it follows also from the "weak" condition that the 
antiquark distributions are not charge symmetric at all points x, but that the integral 
over all x is the same for, say, the up antiquark distribution in the proton and the down 
antiquark distribution in the neutron. 

With the assumptions we have listed, the Gottfried Sum Rule will be equal to 1/3 if 
we have SU(2) flavor symmetry in the proton sea, i.e. if u^ix) = d^ix), or if the integrals 
over X of these distributions are equal. 

One additional point is that the Gottfried "Sum Rule" cannot be obtained from 
current algebra, that is, the GSR cannot be expressed in terms of equal-time commu- 
tators of some observable. The Gottfried Sum Rule is simply a relation which holds in 
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Figure 23: Experimental results for Gottfried sum rule from NMC group, Ref. [Q. Solid 
dots: F| - F^; open dots: - F^{x')]dx' /x'. 



the quark/parton model, with additional assumptions regarding equality of strange and 
charmed quark expectation values in the proton and neutron. This is not the case for 
the Adler sum rule, which can be derived either from current algebra relations or the 
quark/parton model. 

For almost twenty years there have been indications that the GSR seems to be less 
than 1/3. However, interest in the Gottfried Sum Rule has intensified with the very 
accurate measurements in 1991 by the NMC group Q]; the data was re-analyzed in 
1994 [^]. Their results are shown in Fig. The solid circles show F2{x) — F^fx), 



while the open circles plot Sg{x) = /^[FKx) — dx/x. From Eq. |80| we see that the 

difference in F2 structure functions is multiplied by 1/x. This means that measurements 
at low X play a critical role in determining So- Fig. ESI shows that roughly half the 
contribution to the GSR integral comes from the region x < 0.1. The neutron results, 
F2, were inferred from reactions on deuterium. 

Earlier measurements by the SLAG, EMG and BGDMS groups [ |67| , [88| -p0| gave results 
which were lower than 1/3, but these values had errors in Sg of roughly 20%, so that the 
results were within one standard deviation of 1/3. The statistical error in the previous 
experiments was dominated by the lack of data at sufficiently small x. The NMG group 
obtained data for x down to 0.003. The final value from the NMG group was 



5g(0.004 < a; < 0.8) = 0.221 ± 0.008 (stat) ± 0.019 (syst), 



(86) 



The NMG group then fitted the difference in F2 structure functions by a power law and 
extrapolated from x = 0.004 to x = (the extrapolation to a; = 1 produces no measurable 
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contribution to Sg)- Their extrapolated result was 



G 



0.235 ±0.026. 



^7) 



This is more than four standard deviations lower than the "naive" expectation of 1/3. 

If we use the structure functions for electromagnetic interactions, Eqs. |TT| and but 
do not invoke charge symmetry, we obtain 
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A5d{x) + 5u{x) + A5d{x) + 5u{x) 
Integrating this over x and normalizing the valence quarks gives 
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(89) 



There has been much speculation as to the cause of the NMC result. One possibility 
is that the Gottfried Sum Rule is, in fact, 1/3, and that the apparent deviation of the 
GSR from 1/3 is an artifact of the procedure for extrapolating the structure functions 
to X = 0. Martin, Roberts and Stirling |]r2| suggested that one might have Sq = 1/3, 
where the "missing" contribution to Sq comes from very small x values, and that the 
NMC power law extrapolation was in error. 

We discussed in Sect. 4 other theoretical suggestions for the origin of the excess of d^ 
over uP. Recently several groups have compared proton-induced Drell-Yan processes on 
protons and deuterons. Recent results have been obtained for this process by the NA51 
group at CERN fl^, and preliminary results from the E866 group at FNAL |15|. Both 



experiments appear to confirm that d^{x) > vP{x). This was discussed in more detail in 
Section 4. 

Assuming the NMC extrapolation is correct, from Eq. ^ we see that deviation of 
the GSR from 1/3 measures either charge symmetry violation [CSV], or flavor symmetry 
violation [FSV] in the nucleon sea quark distributions (or a combination of the two 
effects). If we assume the validity of charge symmetry, then the NMC measurement 
implies a surprisingly large SU(2) flavor asymmetry in the proton antiquark distributions, 
namely 



dx 



dP{x) 



cs 



0.147 ±0.039 



(90) 



The FSV contribution suggested by the NMC experiment is surprisingly large, as 
it is much larger than can be accommodated by perturbative QCD. Both NLO and 
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NNLO QCD calculations have been carried out, and predict very small FSV effects P7[ . 



Consequently, we need a non-perturbative mechanism to generate ffavor- violating sea 
quark distributions which will reproduce the experimental result. 

The Pauli principle should make it easier to form a dd pair than a uu pair in the 
presence of the valence quarks. In Feynman and Field's early paper on parton distri- 
butions IQ, they assumed an excess of d quarks in the proton on these grounds. A 



promising mechanism for generating additional d quarks in the proton, first recognized 
in Ref. |]SB|, is the "SuUivan Effect" [Q. We discussed this briefly in Sect. 4, and refer 



the interested reader to the review article of Kumano . 

It is important to note that in principle, one could reproduce the NMC results even 
if flavor symmetry is exact] this was pointed out by Ma [jl6[. From Eq. |8^, if we assume 



exact flavor symmetry, but not charge symmetry, in the nucleon sea, then 

dP{x) = uPix) = f{x) 

qP(x) ^ (91) 

then the NMC measurement implies a substantial charge symmetry violation in the nu- 
cleon sea, i.e., 

[ dx \qP{x) - q"ix)] = -0.088 ± 0.023 (92) 
Jo 

It would require a very sizable CSV contribution to reproduce the NMC result. The 
necessary parton sea CSV contribution is more than an order of magnitude larger than 
the theoretical estimate we discussed in Sect. 3.2. Alternatively, the NMC data could 
result from a linear combination of FSV and CSV effects in the nucleon sea. 

Eq. ^ shows that the Gottfried Sum Rule is sensitive to charge symmetry violation 
in the nucleon sea. One can also have charge symmetry violation in the valence quark 
distributions. However, the integral over x of the charge symmetry violating pieces 
must vanish since CSV contributions cannot change the valence quark normalizations. 
In principle, violation of charge symmetry in the valence quark distributions makes no 
contribution to the Gottfried Sum Rule. However, as the valence quark CSV contribution 
vanishes only when integrated over all x, it is possible to obtain a contribution from 
valence quark CSV if data is taken only over a finite range of x. It is important that 
extrapolations over an unmeasured region properly account for these terms. 



6.2 Adler Sum Rule 



The Adler Sum Rule is given by the integral of the F2 structure functions for charged 
current neutrino scattering. The Adler Sum Rule, Sa, can be defined as 



Sa = 



lim 

32^00 Jo 
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Figure 24: Experimental results for Adler sum rule, Eq. |93| from the WA25 group, Ref. 

El- 



= 1 . (93) 

We obtain the result S'a = 1 if we neglect the term |V^dp ^ 1 x 10"'^. The Adler sum 
rule thus requires measuring the F2 structure function for antineutrinos and neutrinos 
on protons, dividing by x (which emphasizes the contribution from very small x), and 
subtracting them. The Adler sum rule then follows from the normalization of the quark 
distributions. As a consequence of the algebra of SU(2) charges, the Adler sum rule has 
no QCD corrections. 

In Fig. |2^ we show the experimental situation regarding the Adler sum rule. The 
experimental data are from the WA25 experiment using the CERN-SPS wide band 
neutrino and antineutrino beams in the BEBC H and D bubble chambers. The experimen- 
tal data are shown for several values of Q^. The average value is Sa = 1-01 ± 0.08 ± 0.18. 
However, as pointed out by Sterman et al. |j2ll], the total vN cross section used by the 
WA25 group is smaller than the presently accepted value If the WA25 value is 

= 1.08 ±0.08 ±0.18. 

The large errors 



readjusted to fit this total cross section their result becomes Sa 

Within the rather large errors, the results are independent of 
arise from the factor 1/a; in the integral, Eq. |9^. This gives a heavy weighting to the data 
at small x. The paucity of data in this region and the relatively large error bars there give 
a large uncertainty in the sum rule value. Because of the difficulties in obtaining sufficient 
neutrino and antineutrino cross sections at small x values for light nuclear targets, it is 
unlikely that there will be new experimental neutrino data in the near future which would 
allow us to test the Adler sum rule. 

Note that the experimental points presented involve neutrino measurements on neu- 
trons (e.g., deuterons) and protons, and not antineutrinos and neutrinos on protons, as 
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given in the definition of tlie Adler sum rule, Eq. ^ Tliis follows from the relation 
j?w+n^^^ q2^ cs pW p^^^ Q'^), which follows if one assumes charge symmetry, as was dis- 
cussed in detail in Section V.A. Thus, the WA25 group does not plot the integral Sa of 
Eq. B^, but instead measures the quantity 



lim 
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6u{x) + 5d{x) 



6d{x) 



(94) 



If charge symmetry is exact, then 5*^ = •S'^- If we assume that the quark normalization 
integral is indeed one, then the experimental results allow us to place an upper limit on 
the integral of the antiquark charge symmetry violating amplitudes at the 20% level. We 
will discuss this more in Sect. 7.2 in connection with the "charge symmetry" sum rule. 



6.3 Gross-Llewellyn-Smith sum rule 



The Gross-Llewellyn Smith [GLS] Sum Rule is derived from the structure functions 
for neutrinos and antineutrinos, 
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which follows from the normalization of the quark valence distributions. We have pre- 
sented the sum rule for an isoscalar target. An identical prediction is obtained for either 



a proton or neutron target 



i.e. xF^"^'^'{x) 



\x 



or xF^ 



(x) + xF._ 



x), in the 



integrand of the sum rule. Neither the Adler nor Gross-Llewellyn Smith sum rules require 
any additional assumptions regarding charge symmetry of quark distributions. The GLS 
sum rule does, however, acquire a QCD correction, which is represented by the term in 
square brackets in Eq. The higher order QCD corrections have been derived by Larin 
and Vermarseren ||93|. They depend on the strong coupling q;s((5^). The terms a and h 
depend on the number of quark flavors (n/) available at a particular value of . The 
quantity IS.HT represents a higher twist contribution ^ . 



As is the case for the Adler and Gottfried sum rules, the Gross-Llewellyn Smith sum 
rule requires that the structure function be divided by x in performing the integral. This 
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Figure 25: Experimental results for Gross-Llewellyn Smith sum rule, Eq. from CCFR 
group, Ref. |^ . 



gives a strong weighting to the small-x region, such that as much as 90% of the sum rule 
comes from the region x < 0.1. The GLS sum rule is the most precisely known of the 
three sum rules we consider. The best (and most recent) value has been obtained by 
the CCFR collaboration |^|, which measured neutrino and antineutrino cross sections 
on iron targets, using the quadrupole triplet beam (QTB) at FNAL. In Fig. ESI we show 
the CCFR measurements, the experimental values of xF^, and their integral, vs. x. They 
obtain cross sections at several values of x and Q^; the final value for the sum rule is given 
for = 3 GeV^. Their reported value for the sum rule is Sgls = 2.50 ± 0.018 ± 0.078. 
The GLS sum rule is therefore known to 3%. 

A theoretical value for the Gross-Llewellyn Smith sum rule requires evaluating the 
QCD corrections. The most recent calculations include next-to-leading order QCD correc- 
tions. They use a QCD scale parameter Aqcd = 213±50 MeV. With this scale parameter 
and NLO QCD corrections, one obtains a theoretical prediction Sgls = 2.63 ±0.04. The 
theoretical prediction is just within two standard deviations of the experimental value. 
In Fig. ^ we show the evolution over time of the GLS sum rule value. 

The errors on the GLS sum rule are now at a level where the value of the strong 
coupling constant ag is a major source of error. The CCFR group may now have data on 
xFs over a wide enough range of that, together with renormalized data from several 
other experiments, they may be able to evaluate the GLS sum rule without extrapolation 
for a large range of values. This raises the hope that one can calculate the Gross- 
Llewellyn Smith sum rule as a function of Q^, and to use the resulting dependence of 
the sum rule to determine as(Q^). 

The CCFR group has recently re-calculated both the GLS sum rule, and the strong 
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CDHS 3.E0 + 0.5 

CHARM 3.58 + 0.43 
COFRR 3.e3 + O.SO 
WA2S 2.70 + 0.40 

CdFR(NBB) 2.70 +0.15 
dC!FH(QTB) 2,50 + O.OS 



Figure 26: Experimental results for Gross-Llewellyn Smith sum rule, and their errors, for 
a series of experiments, in chronological order from top to bottom. 



coupling constant [0]. With data of this quality over a large range, it may be 
possible to use the dependence to put constraints on the strong coupling constant. 
Additional information regarding this procedure can be found in the thesis by Seligman 

The quantities xF^*'^'''[x) + xF^ which form the integrand for the GLS sum 

rule, are obtained by taking the difference between cross sections for neutrinos and an- 
tineutrinos on isoscalar targets. This was discussed in Sect. 2.5. If charge symmetry is 
exact, then the F2 structure functions exactly cancel when we take the difference of neu- 
trino and antineutrino cross sections, and we are left only with the F3 structure functions. 
However, if we allow charge symmetry violation, then insofar as the sum of F3 structure 
functions is defined as the difference between neutrino and antineutrino cross sections on 
an isoscalar target, there are additional contributions to this integrand, i.e. using Eqs. 



30| and l33. 



Ua''''^/dx-da'''^'/dx) = ]-(xFY^''\x,Q^)+xFY'''\x,Q^' 



2G^MnE V ' ' / 2 

= |[2K(x) + <(x)) + 6(s(x)-s(x)) 

- 35u^{x) + 5d^{x)] (96) 

In addition to the valence quark distributions which we obtain in the charge symmetric 
limit, there is an additional contribution of three times the difference between strange 
and antistrange quarks. There is an additional term proportional to the valence quark 
CSV terms. From Sect. 3, the CSV terms are predicted to have opposite signs so their 
contributions should add coherently in Eq. Although we expect charge symmetry 
violating contributions of at most about two percent in the F2 structure functions, the 
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relative contribution to the experimental values of F3 could be somewhat larger. This 
is because in the double differential neutrino cross sections (see, e.g., Eq. |^), the F3 
structure functions are multiplied by a coefficient y — ?/^/2, while F2 has a coefficient 
which is roughly 1. Since the average value of y in these experiments is about 0.2, the 
F2 structure functions will on average be weighted by a factor 5 relative to the F3 terms. 
Thus, naively we expect CSV effects in F2 to be magnified by a factor of about 5 in 
extracting F3. Thus a 2% CSV amplitude in the quark distributions could potentially 
change the extracted value of F3 by as much as 10%. 

Despite the possibility that the F3 structure functions could be modified by quark 
CSV effects, in principle, these CSV amplitudes should have no effect on the GLS sum 
rule. Since the GLS sum rule involves an integral over x, the net contribution to the sum 
rule arising from the CSV effects is 

SS9;IX = t \^d,{x) - 6u,{x)] = (97) 
Jo 

The integral vanishes because the quark valence distributions obey the normalization 
conditions, Eq. |2^. Therefore the integral over x of both 6d^{x) and Su^{x) must be 
exactly zero. In practice, this requires having data over all x, or correctly performing 
an extrapolation over the unmeasured x region, so that the contributions from valence 
quark CSV terms really average to zero. 

7 Flavor Symmetry and Charge Symmetry Sum Rules 



From the previous section, we see that both FSV and CSV terms contribute to the 
Gottfried sum rule, and that what is conventionally called the "Adler sum rule" also 
contains a CSV contribution from sea quarks. If sufficiently accurate experimental data 
can be obtained, it would be useful to derive sum rules which could differentiate between 
charge symmetry and flavor symmetry violation in the nucleon sea. Quantities like the 
electromagnetic interactions, and Drell-Yan processes, which couple to virtual photons 
are proportional to the squares of the quark charges. They will give the up quark (and 
antiquark) distributions a relative weighting four times that for the down quark distri- 
butions. Neutrino deep inelastic structure functions, which couple to the weak isospin, 
allows the possibility of separating these contributions. 

This process of defining sum rules which would differentiate between charge symme- 



try and flavor symmetry violation was originally carried out by Ma |T6[. Our charge 
symmetry and flavor symmetry sum rules are linear combinations of the integrals defined 
by Ma. 
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7.1 Charge symmetry sum rule 



We define a "charge symmetry" sum rule in terms of tlie F2 structure functions for 
charged current neutrino interactions on the neutron and proton, 



Scs = lim / ^ 

2x 



Fr^(a;)-Fr-"(a;) 



W+1 



dx 



dx 



u^{x) + s{x) + rfP(a;) - + s{x) + 



5u{x) + 5d{x) + s{x) — s{x) 



dx 



6u{x) + 6d{x) 



{9t 



In deriving Scs, we assume that the strange and charmed contributions for neutron and 
proton are identicaL We term this the "charge symmetry" sum rule, since by inspection if 
either the strong form or weak form of charge symmetry holds for the nucleon sea quark 
distributions, then Scs will be identically equal to zero, and any deviation from zero 
will be due to violation of charge symmetry. (To reiterate, the "strong form" of charge 
symmetry for sea quarks is the statement that u^{x) = d'^{x) for all x; the "weak form" 
would state that the distributions might not be identical, but that their integrals over x 
are equal. With either form of charge symmetry, the contribution to the integrals in the 
sum rule would vanish). Just as for the Adler sum rule, there are no QCD corrections to 
the charge symmetry sum rule. 

We can relate the charge symmetry sum rule to the sum rules listed in Sect. 6.2. We 
can easily see that 



S 



cs 



S, 



S, 



Sa — Sa where 



lim / dx 

32^00 Jo 



lim / dx 



Fr""(x,g^)-Fr"(x,Q 
2x 



+00 JO 



Fr"(x,g^)-Fr"(x,g^) 

2x 



(99) 



The Adler sum rule Sa requires comparing the structure function F2 for charged current 
weak interactions on protons with antineutrino beams, and with neutrino beams. We dis- 
cussed this at length in Sect. 6.2 (see Eq. |S^). It is required to be one from normalization 
of valence quark distributions, and has no contribution from CSV terms. 

The sum rule Sa requires subtracting the corresponding F2 structure functions for 
neutrinos on neutrons (i.e., deuterons) and protons, respectively; neutrino beams on 
different targets. The "charge symmetry" sum rule requires comparing antineutrinos on 
protons, with neutrinos on neutrons. As an alternative to the charge symmetry sum rule 
one could measure Sa and Sa and compare them. If S'a = 1 and Sa ^ this would give 
clear evidence for charge symmetry violation in the nucleon sea quark distributions. 

One can straightforwardly construct sum rules using different linear combinations 
of F2 structure functions for neutrinos or antineutrinos on protons or neutrons, which 
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contain the same information as the sum rules we have defined here. For example, we 
could define 



Scs = 



lim 



dx 



S'^P = lim 



1 



1 dx 
2x 

6u{x) 

1 dx 
2x 

dx 



X] 



n. 



X\ 



+ 5d{x] 



:,W-p. 



X] 



CS^ 



and 



[X 



Sa + Scs 



6u{x) + 6d{x) 



(100) 



Experimental prospects for accurate measurements of any of these sum rules are poor. 
All these sum rules require precision measurements for both neutrinos and antineutrinos, 
over a wide range of x. Since the sum rules are very sensitive to the small x region, it 
is important to have precise data down to very small x. Furthermore, all the sum rules 
require data on both hydrogen and deuterium targets. As the WA25 group measured 
cross sections from both neutrinos and antineutrinos on protons and deuterium, they 
could in principle construct the charge symmetry sum rule. However, as we have seen 
(viz.. Fig. ^ of Section 6.2), the errors in the Adler sum rule are of the order of 20%, 
so the charge symmetry sum rule would be consistent with zero at the 20% level. We 
discussed in the previous section the difficulties in obtaining precise neutrino data on 
light targets, over a wide range of x. 

We can compare our "charge symmetry" sum rule with the one proposed by Ma [16 
He defined the following sum rule (Eq. [14] of Ref. 

"1 dx 



m) 



~,{Ma) 
^CS 



2x 



X 



+ Fr'(x) 



W+r 



(x) + R 



W-ni 



\X 



(101) 



We see that if the charge symmetry sum rule can be written as S'qI^"'^ = Scs + Scs', and 



since we have shown that Scs 
information. 



Scs, all three of these sum rules give precisely the same 



7.2 Flavor symmetry sum rule 



We can define a "fiavor symmetry" sum rule by comparing the structure functions 
from charged current neutrino interactions on protons and neutrons, i.e. 

dx r 



Sfs = 



xFr''-ix)-xF,^'^ix)\ 
1 

dx[— {U^{x) + c(x) — d^{x) — s{x)) 



+ rf"(x) + s{x) - m"(x) - c{x)] 



dx 



- dlix) + u^ix) - m"(x) - d^{x) + c/"(x) 

\ asm 



dx 



uP{x) - dP{x) - u"(x) + rf"(x) 



TT 



(102) 
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If flavor symmetry holds for the nucleon sea, i.e. if u^{x) = dP{x) and u^{x) = (or 
the weak condition that the integrals over x of these distributions are equal), then Sps 
will be equal to one. So any deviation of this sum rule from one signifies flavor symmetry 
violation in the nucleon sea. There is a QCD correction to this sum rule. It is the same 



as for the Gross-Llewellyn Smith sum rule (see Section 6.3). In Eq. 102 we include the 
lowest-order QCD correction. 

Experiments with muons suggest the magnitude of flavor symmetry breaking which 
we might expect. As we discussed in Sect. 3.2, theoretical calculations predict that sea 
quark CSV effects will be very small. In that case the antiquark contributions to Eq. 
102| will be identical to those measured by the NMC group [Q. For = 10 GeV^, we 



then expect the integral over the antiquark distributions to give a 30% effect. However, 
there will be significant experimental difficulties obtaining accurate data for the sum 
rule. The factor 1/x in the integrand requires precise neutrino data on both protons 
and neutrons (i.e., deuterium), at very small x. Even small differences between the 
structure functions on protons and deuterium become magnified when one integrates the 
difference between them. 

Our flavor symmetry sum rule was obtained by combining structure functions for 
neutrinos on protons and neutrons. We could define an analogous function by utilizing 
antineutrino structure functions on protons and neutrons. 



Sfs = 



1 dx r 



2x 
1 

dx 





xF^~P{x) -xF^'^'i 



X 



<(x) - <(x) + uP{x) - m"(x) - cF(x) + d"(x) 



1+ dx uP(x) - dP(x) - m"(x) + rf"(x) 

Jo L 



. (103) 



Comparing Eq. p.03| with Eq. p.02| we see that Sps = Sps] consequently both of these sum 
rules give exactly the same information. 



S' = 



Ma defines an additional sum rule, S' (Eq. [23] of Ref. |1T6[|), through 

/■I dx 
lo X 

asm 



xF^''P(x) - xFf "^(x) - xFf ^"(x) + xFf 



X] 



-4 - 4 / dx 
Jo 



uP{x) - dP{x) - m"(x) + rf"(x) 



1 - 



IT 



(104) 



This is related to the two previously defined sum rules by S' = ~2{Sfs + Sfs), and 
gives the same information as either of those. It requires knowledge of both neutrino and 
antineutrino parity violating structure functions on both the proton and neutron. 



8 Conclusions 
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We have reviewed the vahdity of charge symmetry for parton distributions. We calcu- 
lated CSV contributions for both valence and sea quarks. The "majority" valence quark 
distribution, and the sea quark distribution, are predicted to obey charge symmetry to 
at least 1% for all values of x. The "minority" valence quark distribution, however, is 
predicted to show CSV effects of between 3-7% at large Bjorken x. This prediction ap- 
pears to be robust, as all calculations of quark CSV predict effects of this magnitude, 
and reasonably model-independent estimates confirm these results. 

This violation of charge symmetry is large. If this turns out to be correct, we should 
probably re-evaluate all phenomenological parton distributions, introduce some explicit 
charge symmetry violation and re-fit existing cross sections. Towards this end, we have 
redefined a number of observables, using a formalism which does not assume explicit 
parton charge symmetry. One difficulty in searching for CSV effects at present is that 
all phenomenological parton distributions assume charge symmetry at the outset, so any 
existing CSV effects have been absorbed into the current parton distributions. This makes 
it difficult to search for experimental violation of charge symmetry. Of course, if CSV 
effects were extremely large, we would already have seen this in existing experimental 
data. 

If charge symmetry is indeed violated at the predicted level for the minority valence 
quarks, then one should be able to measure such effects. First, we reviewed the status of 
current tests of parton charge symmetry. There are essentially two such tests. The first 
involves the comparison of charge-current cross sections induced by neutrinos to those 
from antineutrinos. This comparison can detect CSV in valence quark distributions; the 
test is only valid at large x. Existing data gives only weak upper limits on valence quark 
CSV. The second test of parton CSV is the so-called "charge ratio," the ratio of F2 
structure functions in DIS processes induced by leptons, to the F2 structure functions in 
charge- changing weak processes from neutrino beams (both of these need to be measured 
on isoscalar targets). Despite the fact that this comparison requires a large number of 
corrections, it appears that most of these corrections are under control. 

Recent data, muon induced DIS on deuterium from the NMC group and neutrino 
reactions on iron from the CCFR group at FNAL ||6^, allow us to place rather tight 
constraints on parton CSV amplitudes. Such relations should hold at all values of x. 
These experiments allow us to set limits of roughly 10% on charge symmetry violating 
contributions to parton distributions, in the region 0.2 < x < 0.4. However, in the region 
X < 0.2, there appears to be a discrepancy between the electromagnetic and weak F2 
structure functions. 

In this review, we have suggested several experiments specifically designed either to 
detect CSV in parton distributions, or to set more stringent upper limits on parton CSV. 
Probably the most sensitive test would be a comparison of the F2 structure functions 
for charge- changing weak processes induced by electrons and positrons on an isoscalar 
target. At sufficiently high energies (such as are accessible at HERA), comparison of 
— D and e~ —D processes should be a powerful and relatively clean test of parton CSV. 
Other potentially useful reactions which are quite sensitive to parton CSV are Drell-Yan 
processes for charged pions on isoscalar nuclear targets, or semi-inclusive leptoproduction 
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of pions on isoscalar targets. We made estimates of the magnitude of CSV effects expected 
in these reactions. 

We also discussed the contributions of CSV effects in various sum rules. We showed 
that CSV contributions to the nucleon sea have no effect on the Gross-Llewellyn Smith 
sum rule, but in principle they affect both the Adler and Gottfried sum rules. In fact, the 
existing experimental "test" of the Adler sum rule contains a contribution from CSV in 
sea quark distributions, and can be used to place an upper limit on parton sea quark CSV. 
We introduced two new sum rules, a "charge symmetry" sum rule which is zero if charge 
symmetry is exact, and a "flavor symmetry" sum rule which is one if flavor symmetry 
is exact in the nucleon sea. Both sum rules require measuring the structure functions 
for neutrinos and/or antineutrinos on isoscalar targets, with particular attention to the 
small-x region. Limits on CSV from present experiments are in the neighborhood of 10% 
of the average parton distributions. 
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